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In the baryogenesis via leptogenesis scenario the self-energy contribution to the CP-violating pa- 
rameter plays a very important role. Here, we calculate it in a simple toy model of leptogenesis using 
the Schwinger-Keldysh/Kadanoff-Baym formalism as starting point. We show that the formalism 
is free of the double-counting problem typical for the canonical Boltzmann approach. Within the 
toy model, medium effects increase the CP-violating parameter. In contrast to results obtained 
earlier in the framework of thermal field theory, the medium corrections are linear in the particle 
number densities. In the resonant regime quantum corrections lead to modified expressions for the 
CP-violating parameter and for the decay width. Most notably, in the maximal resonant regime the 
Boltzmann picture breaks down and an analysis in the full Kadanoff-Baym formalism is required. 
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I. INTRODUCTION 

From the theoretical point of view the baryogenesis 
via leptogenesis scenario [1] is a very attractive explana- 
tion of the observed baryon asymmetry of the Universe. 
One of its key ingredients is heavy Majorana neutrinos, 
whose CP-violating decays are responsible for the gen- 
eration of a lepton asymmetry, i.e. for leptogenesis. The 
CP-violating parameter receives contributions from the 
vertex [1, 2] and self-energy [3-6] diagrams. If the masses 
of the heavy neutrinos are strongly hierarchical, then the 
two contributions are comparable. However, if the mass 
spectrum is quasidegenerate then the self-energy contri- 
bution is resonantly enhanced and becomes considerably 
larger than the one from the vertex diagram. The reso- 
nant enhancement of the CP-violating parameter allows 
one to bring the scale of leptogenesis down to ~ 1 TeV 
[6, 7]. This scenario is very interesting from the experi- 
mental point of view since it is potentially accessible in 
accelerator experiments [8]. 

In state-of-the-art calculations the self-energy CP- 
violating parameter is evaluated in vacuum and then used 
to calculate the asymmetry generated by the decays of 
the heavy neutrinos in the hot and dense medium. In this 
approximation the possibly important medium effects are 
neglected from the very beginning. One can take them 
into account systematically by using the Schwinger- 
Keldysh/KadanofT-Baym formalism [9] or approximate 
self-consistent equations derived from the Kadanoff- 
Baym equations (see [2, 10-22] for related work). 

The Kadanoff-Baym formalism provides a powerful 
framework for studying nonequilibrium processes within 
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quantum field theory. However, it is technically consid- 
erably more involved than the canonical Boltzmann ap- 
proach. For this reason, we apply it here to a simple toy 
model which we have already used in [2] to investigate 
the vertex contribution. The Lagrangian contains one 
complex and two real scalar fields: 
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C = -d^id^i - -Mt^i + d^b - m'bb 
- ^(bb) 2 - |Vi66 - ^ipibb + C rest ,i = l,2. 



(1) 



where b denotes the complex conjugate of b. Here and 
in the following we assume summation over repeating 
indices, unless otherwise specified. Despite its simplic- 
ity, the model incorporates all features relevant for lep- 
togenesis. The real scalar fields imitate the (two light- 
est) heavy right-handed neutrinos, whereas the complex 
scalar field models the baryons. The U(l) symmetry, 
which we use to define "baryon" number, is explicitly bro- 
ken by the presence of the last two terms, just as the B— L 
symmetry is explicitly broken by Majorana mass terms 
in phenomenological models. Thus the first Sakharov 
condition [23] is fulfilled. The couplings <?.; model the 
complex Yukawa couplings of the right-handed neutri- 
nos to leptons and the Higgs. By rephasing the complex 
scalar field at least one of the couplings gi can be made 
real. If arg(gi) ^ arg(^2) the other one remains complex 
and there is CP-violation, as is required by the second 
Sakharov condition. In vacuum the self-energy contribu- 
tion to the CP-violating parameter is given by: 
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where Tj is the decay width of the heavy scalar ipj. 
Note that the expression (2) has been obtained using 
the formalism developed in [6] . The formalism employed 
in [24] leads to a slightly different expression for the 
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CP-violating parameter (see Appendix A for more de- 
tails). The required deviation from thermal equilibrium 
is caused by the rapid expansion of the Universe. Thus 
the third Sakharov condition is fulfilled as well. Finally, 
the quartic self-interaction term in (1) plays the role of 
the Yukawa and gauge interactions in established mod- 
els - it brings the "baryons" to equilibrium. The renor- 
malizability of the theory requires the presence of some 
additional terms, which are accounted for by £ re st- 

Here, we study the self-energy contribution to the CP- 
violating parameter in the hierarchical and quasidegener- 
ate cases using the Schwinger-Keldysh/Kadanoff-Baym 
formalism as starting point. To make the discussion less 
technical we give the details of the calculation in the 
appendixes, whereas in the main body of the paper we 
sketch the derivation and present the results. 

(i) As we argue in Sec. Ill, the formalism is free of 
the double-counting problem typical for the canon- 
ical Boltzmann approach. In other words the struc- 
ture of the equations automatically ensures that the 
asymmetry vanishes in thermal equilibrium and no 
need for the real intermediate state (RIS) subtrac- 
tion arises. This property has already been ob- 
served in case of the vertex contribution. 

(ii) The medium corrections to the CP-violating pa- 
rameter are only linear in the particle number den- 
sities. That is, our result differs from that obtained 
previously in the framework of equilibrium thermal 
field theory by replacing the zero temperature prop- 
agators with finite temperature propagators in the 
matrix elements used in the Boltzmann equation. 

(iii) For scalars the medium effects always increase the 
CP-violating parameter, which in turn leads to an 
enhancement of the generated asymmetry. 

(iv) The canonical expression for the CP-violating pa- 
rameter is only applicable in the hierarchical case 
even though it docs not diverge in the limit of equal 
masses. For quasidegenerate masses one has to take 
into account quantum corrections to the effective 
masses and decay widths of the heavy particles in 
medium, which leads to a modified expression for 
the CP-violating parameter. 

(v) In the resonant regime quantum corrections also 
lead to an enhancement of the total in-medium de- 
cay widths. This results in a faster decay of the 
heavy particles and can increase the importance of 
the washout process. 

(vi) In the maximal resonant regime the Boltzmann pic- 
ture breaks down and we argue that an analysis in 
the full Kadanoff-Baym formalism is required. 

In Sec. IV we present numerical solutions of the quantum- 
corrected Boltzmann equations, and discuss the quanti- 
tative impact of medium effects on the final asymmetry 
within the toy model. Finally, in Sec. V, we summarize 
the results and present our conclusions. 



II. KINETIC EQUATION 

The canonical approach to the calculation of the asym- 
metry generated at the epoch of leptogenesis is based on 
the use of Boltzmann equations. In the expanding Uni- 
verse the Boltzmann equation for the toy-baryon distri- 
bution function can be written in the form [25] 

p a V a f b = i[S < (l + / b )-/ b S > ] , (3) 

where all functions are evaluated at the same point (X, p) 
of the phase-space, T> a is the covariant derivative and the 
quantities Eg correspond to the gain and loss terms. For 
decays into a pair of toy-baryons they are given by 

E< (X, p) = - J dn 3 q dnl (2ir) 4 6(k -q-p) 

x \M\l^ bb USX,k)[l + f b (X iq )}, (4a) 

e> (x, p) = - J dn^dnl (2ir)H(k -q-p) 

x \M\ 2 bb ^J b (X,q)[l + f^(X,k)}, (4b) 

where dU^ = d 3 p/[(27r) 3 2E] is the invariant momentum- 
space volume element. The analogous equations for the 
toy-antibaryon distribution function f b and the corre- 
sponding quantities E^ can be obtained from (3) and 
(4) by replacing the subscript b with 6. 

At tree level, see Fig. 1 (a), the decay is CP-conserving, 
that is \M\%^ bb = \M\ 2 bb ^ = \M \ 2 . There are 
two distinct contributions to the CP-violating param- 
eter, the vertex [1, 2] and the self-energy [3, 5] ones. 
The leading-order self-energy contribution is generated 
by the one-loop self-energy diagram depicted in Fig. 1 
(b). Only the diagram where the initial (ifii) and inter- 




(a) (b) 
FIG. 1: Tree-level and one-loop self-energy diagrams of the 
decay process tpi —> bb. 

mediate (ipj) toy-Majoranas are different (i ^ j) con- 
tribute to the CP-violating parameter. Summing the 
tree-level and the "off-diagonal" one-loop self-energy am- 
plitudes one obtains \M\ 2 h .^ bb = |.Mo| 2 (l + £ i) and 
|A^lbfc-s.Vi = |A^o| 2 (l — ei). In the corresponding expres- 
sions for the "antibaryons" the signs in front of ti are 
reversed. If these amplitudes are now substituted into 
E^ and E^ one finds that a nonzero asymmetry is gen- 
erated even in thermal equilibrium, which is inconsistent 
with the CPT symmetry. This is a manifestation of the 
so-called double- counting problem typical for the canon- 
ical Boltzmann formalism. Let us briefly explain what 
that means. An inverse decay bb — ¥ %pi immediately fol- 
lowed by a decay ipi — ¥ bb is equivalent to the two-body 
scattering process bb — ¥ ?pi — ¥ bb where the intermediate 
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toy-Majorana is on the mass shell. That is, the same con- 
tribution is taken into account twice: once when the de- 
cay and inverse decay processes are considered and once 
when the two-body scattering processes are considered. 
This problem is usually solved by subtracting the contri- 
bution of the on-shell intermediate state to the scattering 
amplitude [5, 6]. Roughly speaking, after the subtraction 
the amplitude of the inverse decay process bb —¥ ipi also 
becomes proportional to (1 + ej), which ensures that in 
thermal equilibrium no asymmetry is produced due to 
detailed balance. 

In the canonical bottom-up approach, which has been 
outlined above, one uses elements of the S'-matrix (in- 
out formalism) to calculate the functions S^. In con- 
trast to that, in the top-down approach based on the 
Schwinger-Keldysh/Kadanoff-Baym [26-39] formalism, 
the functions can be identified with self-energies and 
are derived, using nonequilibrium held theory techniques 
(see e.g. [38, 40]), from the two-particle- irreducible (2PI) 
effective action 1 [42] formulated on the closed real-time 
path (in-in formalism). A two- loop contribution to the 




(a) (5) (c) 



FIG. 2: Two-loop contribution to the 2PI effective action and 
the corresponding contributions to the self-energies of the real 
and complex fields. 

effective action and the corresponding contribution to the 
self-energy of the "baryons" are presented in Fig. 2 (a) 
and Fig. 2 (c) respectively. Wigner-transforming the self- 
energy we obtain [2] 

(x, P ) = - J diikdn^TT^s 9 (k - q - P ) 

x g*g 3 Gl(X,k)D $ (X,q), (5a) 
£%(X,p) = - f dIL k dTl q (2Tr) 4 5 9 (k - q -p) 

X g i9 ^G^(X,k)D^(X,q), (5b) 

where dH p = d 4 p/(27r) 4 , S 9 (p) = y r ^g x S(p) is the co- 
variant generalization of the Dirac ^-function and G> 
and are the Wightman propagators of the real and 
complex fields respectively. In the Boltzmann-limit the 



1 Note that there exists a straightforward generalization to a hier- 
archy of so-called nPI effective actions [41]. fn general, for sys- 
tems far from equilibrium and for strong couplings the 2-body 
decay is best described using 3PI. However, for leptogenesis, 2PI 
is sufficient since the CP-violation is dominated by the leading 
loop diagrams. 



latter are related to the distribution function by 

D < =f b D p , D > = (l + f b )D p , (6) 

where in the quasiparticle approximation (see [2] for more 
details) the spectral function reads 

D p (X lP ) = 27rsign(p ) 5 (g^p^pu - m 2 ) . (7) 

In the same approximation the Wightman propagators 
of the toy-antibaryons, D^, can be obtained from (6) by 
replacing f b with f b . 

If the real scalar fields would not mix, i.e. if the off- 
diagonal components of the Wightman propagator were 
equal to zero, we could also write analogous expressions 
for the diagonal components: 

Gii r /~tii /i , r \ /~fii / q\ 

< — H^ ! (j > ~ \ V + Hil^p ' (°) 

with the spectral function given by 

Gf(X,p) = 27rsign(p ) S (g^p, P „ ~ M?) . (9) 

Because of the presence of the Dirac (5-function in (7) 
and (9) the integrations over k° and q° could then be 
performed trivially and we would recover (4) in the CP- 
conserving limit. 

A nonzero toy-baryon asymmetry can be generated 
only if ^ E^. Comparing (5a) and (5b) we see 
that if initially the system is symmetric, i.e. f b = f b , 
this amounts to the requirement that the Hermitian ma- 
trix G^s has complex off-diagonal components. If the 
off-diagonal components peak on the mass shell of the 
quasiparticle species, i.e. if they can be represented in 
the form 

Gl= Si G%+e*G» (i^j), (10) 

then the generation of the asymmetry can be analyzed 
in terms of CP-violating parameters. Substituting (10) 
into (5) we find 

ei = -2Im(5j/5i)Imei . (11) 

To calculate the decomposition coefficients £j we use 
the nonequilibrium formulation of the Schwinger-Dyson 
equation, which is discussed in the following section. 

III. CP- VIOLATING PARAMETER 

Just like in the canonical analysis [5, 6] the starting 
point of our analysis is the Schwinger-Dyson equation 

[G^ 1 ] ij (x, y) = ^-^{x^-W^y), (12) 

where G % i is the full dressed propagator of the "heavy 
neutrinos" , ^ y is the diagonal propagator of the free 
fields and II^ is the self-energy. In the two-loop approx- 
imation, see Fig. 2, the self-energy is given by 

W\x,y) = -\g l9 *D 2 {x,y) - \g*g J D 2 {x,y) , (13) 
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where D(x, y) is the full propagator of the complex scalar 
field and D(x,y) = D(y,x). Note that the arguments of 
the two-point functions and self-energy are defined on the 
positive and negative branches of the Schwinger-Keldysh 
closed real-time contour C [26-31] shown in Fig. 14. 

If wc multiply (12) by G from the right, decompose 
the full propagator and the self-energy into statistical 
and spectral components and integrate over the contour, 
we obtain a system of so-called Kadanoff-Baym equa- 
tions for the statistical propagator and spectral func- 
tion (see Appendix B). Within leptogencsis, the light 
particle species are usually assumed to be very close to 
kinetic equilibrium. Translated to the toy model, this 
means that the toy-baryons described by the propaga- 
tors D(x, y) and D(x, y) are close to equilibrium. In this 
case, one can approximate the full Kadanoff-Baym equa- 
tions by quantum kinetic equations, and find approxi- 
mate analytic solutions of these equations. Using these 
solutions, it is possible to explicitly obtain the decompo- 
sition coefficients £j, which then yield the corresponding 
CP- violating parameters Ci as described above. This ap- 
proach is pursued in Appendix B. In this section we will 
use an alternative approach, which essentially relies on 
the same assumptions, but is somewhat more elegant. 

In the following, we use a compact matrix notation 
for Eq. (12), where we denote the matrices by a hat, 
e.g. G = (G 4 - 7 ). Let us split the self-energy matrix n 
into the diagonal 77 and off-diagonal 77 components and 
introduce a diagonal propagator Q defined by the equa- 
tion 

g- 1 (x,y)=^- 1 (x,y)-n(x,y). (14) 

Subtracting (14) from (12) we find 

G- 1 (x,y)=g- 1 (x,y)-ff(x,y). (15) 

Multiplying (15) by G from the left, by Q from the right 
and integrating over the contour C we obtain a formal 
solution for the full nonequilibrium propagator: 2 

G(x,y) = 6{x,y) 

+ JJ S> A u 3) 4 v G(x, u)ff (u, v)G(v, y) . (16) 
cc 

The invariant volume element, 2# A u = y/—gd 4 u , where 
g = detg^, ensures that (16) can be applied to 
the analysis of out-of-equilibrium dynamics not only in 
Minkowski, but also in a general curved space-time [19]. 
Using the decomposition 

G(x,y) = G F (x,y) - § sign^ - y°)G p (x,y) , (17) 



2 Here, we implicitly assume that the correlators arc diagonal at 
the initial time of the closed time path evolution. Since we con- 
sider the kinetic limit later on, for which the initial time is for- 
mally sent to negative infinity, the results arc not affected. 



and an analogous relation for the self-energy we can split 
(16) into spectral and statistical components. These are 
related to the Wightman propagators by 

G^(x,y) = G F (x,y)T^G p (x,y). (18) 

Because of the sign c function in the decomposition of the 
two-point function and self-energy the integrals over the 
closed-time-path contour in (16) reduce to integrals over 
parts of the uv plane; see Appendix C. Introducing the 
retarded and advanced propagators, 

G R (x,y) = e(x -y°)G p (x,y), (19a) 
G A (x,y) = -6(y - x°)G p (x,y) , (19b) 

and also the retarded and advanced self-energies we can 
represent the resulting expressions as integrals over the 
whole uv plane. Finally, building the linear combinations 
(18) we find for the Wightman propagators 

(x, y) = g^(x, y) - ff&u&v 9(u°) 9(v°) 
x [G R (x,u)n^(u,v)g A (v,y) 
+G^(x, u)n A (u, v)Q A (v, y) 
+G R (x,u)n R (u,v)g^(v,y)]. (20) 

Using (18) and (19) we can also derive formal solutions 
for the retarded and advanced propagators from (20). 
They read 

G R{A) (x, y) = g R(A) (x, y) - ff2>*u$> 4 v 9(u°) 9(v°) 

x G R{A) (x, u)ff R{A) (u, v)g R(A) (v, y) . (21) 

In thermal equilibrium the two-point functions and self- 
energies in (20) and (21) depend only on the relative co- 
ordinate, s = x — y, and are independent of the center co- 
ordinate, X = \ {x + y), i.e. are translationally invariant. 
As discussed above, we expect that the deviations from 
equilibrium are moderate. In this case, one can perform 
a gradient expansion of the two-point functions and the 
self-energies in the vicinity of X keeping only the leading 
terms. Performing the Wigner transformation, i.e. the 
Fourier transformation with respect to the relative coor- 
dinate (see Appendix B) wc trade the relative coordinate 
s for a coordinate p in momentum space. Effectively, 
the Wigner transformation replaces the coordinate-space 
arguments (x, y) of each two-point function in (20) and 
(21) by the phase-space coordinates (X,p) and "removes" 
the double integration; see Appendix C. Combining the 
Wigner transforms of (20) and (21) we find for the full 
Wightman propagators 

^ [i gA] - g R ff^g A ] [i ff A g A ] 
* det[i-g R n' R ]det[i-n' A g A ] 

where all the functions are evaluated at the same point 
{X,p) of the phase-space. The term proportional to 77^ 
describes scattering [37-39, 43-48] and is irrelevant for 
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us at the moment. Contracting the decay term with the 
product of the couplings in (5a) and using that, by con- 
struction, the propagators Q are diagonal and 77' off- 
diagonal, we find 



9i9j 



det[i-g R n' R ] det[i-n A g A ] ' 



(23) 



where 



| 5 ,| 2 A* ee \gtfTl^G%G% - 2Re[ 9i g*g%Il%} . (24) 

In (24) and all the equations below we implicitly assume 
summation over all j ^ i. 

A very important feature of (23) is that the loop cor- 
rections A^ are the same for both the gain and loss terms 
(i.e. for the > and < components), respectively. To ob- 
tain an equivalent result in the canonical approach one 
needs to apply the real intermediate state subtraction 
procedure [5, 6]. This means that, here, the structure of 
the equations automatically ensures that no asymmetry 
is generated in thermal equilibrium. Stated differently, 
the formalism is free of the double- counting problem and 
no need for RIS subtraction arises. This conclusion is in 
accordance with the corresponding result for the vertex 
contribution to the CP-violating parameter [2]. 

Contracting the decay term with the product of the 
couplings in (5b) we obtain an expression similar to (23) 
but now with 



\9A 2 rtMe>M - 2Re[g*g 3 g%n%\ . (25) 



From Eq. (5), we see that the rates for the decays ipi — > bb 
and ipi — > bb differ from each other only if ^ A|. The 
corresponding CP- violating parameter is given by 



£i = 



[1 + A*]-[1 + A|] 
[1 + A*] + [1 + Ai] 



(26) 



The first terms in (24) and (25) are equal and cancel out 
in the difference of A^ and A|. Therefore, in agreement 
with (11), the CP-violating parameter can be expressed 
as 



£i = 



-21m( gj -/ gi )lm(gffll£) 

1 + H A t + A l) 



(27) 



Equation (23) also provides us with an expression for 
the total in-medium decay widths of the heavy particles: 



-^rned 



= r, x 



1 + ±(A«+A^ 



det [i - g R n R ] det [i - n' A g A ] 



(28) 



where Ti — \gi\ 2 / (16irMi) are the corresponding tree- 
level decay widths in vacuum. 

The CP-violating parameter (27) is a function of the 
space-time coordinate X and four-momentum p. Note 
that in general p does not need to be on-shell. The on- 
shell condition is determined by the diagonal components 



of the full propagator. The CP-violating parameter car- 
ries information about the off-diagonal components, com- 
pare Eq. (10). Therefore, the conventional, "on-shell" , in- 
terpretation of (27) is only applicable if the off-diagonal 
components of the full propagator peak at the same val- 
ues of the momentum as the diagonal ones. Furthermore, 
if we want to use the Boltzmann equations for Q to calcu- 
late the asymmetry, then the mass spectra of the diagonal 
(g) and full (G) propagators must be sufficiently close. 

The value of the CP-violating parameter depends on 
the masses Mi of the heavy species and the coupling con- 
stants gi or, alternatively, the decay widths T^. It is use- 
ful to discriminate between three cases, namely 

M\ < M|, M 2 > MjTj (strongly hierarchical case) , 

(resonant case) , (29) 
(maximal resonant case) , 



Ml »|AM2.|»M J T j 
Mf »|AM 2 .|<M j r j 



M? — M?. As we will argue in the fol- 



where AM?- 

lowing, the generation of the asymmetry can be studied 
approximately in the Boltzmann picture using the above 
CP-violating parameter in the strongly hierarchical and 
resonant cases, but not in the maximal resonant case. 



A. Strongly hierarchical case 

In this subsection, we derive the CP-violating param- 
eter in the strongly hierarchical limit, Mi -C M2. As 
we shall see later, in the hierarchical case the denomi- 
nator in Eq. (27) is close to unity. The same is true for 
the denominator of (28). Therefore the in-medium de- 
cay widths r™ ed coincide with the tree-level ones. For 
the CP-violating parameter we obtain 



-2Im( 5j / 5i )lm(^n£) 



(30) 



It only remains to calculate the last term in (30). Let 
us first analyze the self-energy YIr. In a "baryonically" 
symmetric configuration the particles and antiparticles 
are interchangeable, so that D(x,y) = D(x,y) (see Ap- 
pendix F for a proof). In the following, we will use that 
the baryon asymmetry is small, i.e. D(x,y) ~ D(x,y). 
In this case the two terms of (13) combine to a sym- 
metric matrix which is proportional to (g* gj + gig*)- To 

obtain Hr(X,p) we decompose this matrix into a statis- 
tical and spectral part, multiply the spectral part with 
9(x° — y°) and perform the Wigncr transformation. The 
result reads 

tt%(X,p) = - (g*g 3 + g^Jdllkdll^nfSip - k - q) 

(31) 



xD$>(X,k)D%>(X,q), 

where the superscript 's' refers to the "baryonically" sym- 
metric configuration and Dp is the statistical propagator 
of the complex field. In the quasiparticle approximation 
Dp = (i + fb)D p . The retarded and advanced propa- 
gator of the complex scalar field can be represented in 
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the form Dr(a) = Dh ± § D p , where in the quasiparticle 
approximation D p is given by (7) and Dh 



v 

■> '■> ; sec 



Eq. (G7). Consequently, in the symmetric configuration 
the retarded self-energy is given by a linear combination 
of two real- valued symmetric matrices: 

iTi = n« + jn« = + .g 4 .g*)(L, + • L„), (32) 

where, to shorten the notation, we have introduced 

L h{p) (X, p) = Wn Jdn k dn q (2tt) 4 % - k - q) 

x D { ;\x,k)D h s l p) {X,q) . (33) 

For a homogeneous and isotropic system the one-particle 
distribution functions depend only on the Lorentz- 
invariant product ku, where k is the particles' momentum 
and u is the (constant) four-velocity of the medium with 
respect to the chosen frame of reference. Using this prop- 
erty, in Appendix G we evaluate L p on the mass shell of 
the z'th quasiparticle: 

L p (X,p) = rj f [1 + f b (E p ) + f- b (E p )] , (34) 



.,, = \[Ei + r\p\ cos6]. Here, E t = (M? + p 2 )" 



where E. t 

and p are the components of the four-momentum p of 
the heavy scalar in the rest-frame of the medium, and 
r = (1 — 4m 2 /M 2 ) 2 . For massless toy-baryons in vacuum 
L v p ac \ m =o = 1. Taking a nonzero toy-baryon mass m into 
account yields L v p ac — r. In medium (symmetric case), 
we find using Eq. (34) 



L P (X,p) 

Jjuac 



=i 



2 r^ma 

T IpI JE mtn 



= /2 

fb(E)dE, 

/2 



(35) 



where E max = Ei + r\p\ and E min = E { - r|p| are the 
largest and smallest kinematically allowed energies of the 
light scalars produced in the decay ipi — > bb. 

As we show in Appendix B, analogously to (32) the 
Wigner transforms of the retarded and advanced propa- 
gators can be represented as linear combinations of two 
Hcrmitian matrices. Applied to the matrix Q this implies 
that it splits into two real- valued diagonal matrices: 



Q33 _ Q33 i iQjj 
y R(A) ~ y h ^ 2 y p 



(36) 



In the hierarchical case we can neglect the finite decay 
widths of the quasiparticle species and approximate Q 3 . 3 
by (9). It diverges on the mass shell of the j'th quasi- 
particle and is zero everywhere else. Because of the pres- 
ence of G^ in (23), the CP- violating parameter (30) must 
be evaluated on the mass shell of the z'th quasiparticle. 
Since the product of G*| and Q 33 vanishes in the quasipar- 
ticle approximation, we conclude that only the Gh term 
contributes. It is given by [19] 



In vacuum the tree-level decay width of the heavy species 
is given by 1^ = \gi\ 2 / (WwMi). Therefore, we can rewrite 
the diagonal components of the spectral self-energy in the 
form W p l = -2MiTiL p . 

In the hierarchical case | M 2 - M 2 1 > | Iljj - 11^ | . Sub- 
stituting (32) and (36) into (30), and evaluating the mo- 
mentum on-shell, p 2 — M 2 , we then obtain for the CP- 
violating parameter 3 



e i = -Tm( 9j /g i )g> h >ILi i 



\9j\ 
16tt 



dm 



fJO-J) 



M 2 



Mf 



gfgJiMf-Mf^ + iMjTj^lA 



(38) 
■L . 



For massless toy-baryons in vacuum L p = 1 and we re- 
cover the classical result (2). In the strongly hierarchical 
limit, one can neglect the decay width Tj in the denom- 
inator of (38). In this case, 



e i /eV a ° = L p {X,p)/L': 



(39) 



This result is identical to the one for the vertex contri- 
bution [2] . The CP-violating parameter is given by a sum 
of vacuum and medium contributions. From Eq. (35), we 
find that the medium contributions are proportional to 
the one-particle distribution function, which is positive. 
Hence, for the scalar toy model and the strongly hierar- 
chical limit considered here, the self-energy contribution 
to the CP-violating parameter is always enhanced by the 
medium effects. 

Note that the medium contribution depends only on 
the distribution function of the toy-baryons and is inde- 
pendent of that of the toy-Majoranas. Since we expect 
the light scalars to be close to kinetic equilibrium at all 
times, it is instructive to estimate the size of the correc- 
tions in thermal equilibrium. In the hierarchical limit the 
asymmetry is predominantly generated by the decay of 
the lighter toy-Majorana. Inserting a Bose-Einstein dis- 
tribution (BE) or a Maxwell-Boltzmann (MB) distribu- 
tion function we obtain for the ratio of the corresponding 
CP- violating parameter and its vacuum value: 



ei(|p|) 



cxp ( 


E, 


2T 


cxp ( 


E 


2T 






ax -2/i 


— e 




2T 



BE , 



MB . 



The temperature- and momentum dependence of the 
medium correction in the range of typical momenta |p| ~ 
T is shown in the shaded areas in Fig. 3 for the BE and 
MB cases, respectively. We also show the CP- violating 
parameter (ei) obtained by averaging Eq. (38) over the 
momentum |p|. As expected, (ei) ~ ei (|p| <~ T). 



p- 



Mf - IE* 



(p 2 - M 2 



n u )2 



n«)s 



3 Note that, by setting p 2 = M 2 , we neglect possible off-shell 

(37) 

effects (sec Appendix A). This is a good approximation if |M; — 
Mj \ 2> Fj, which is the case in the strongly hierarchical case. 
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FIG. 3: Effective self-energy CP-violating parameter ei in 
medium obtained from the Kadanoff-Baym formalism. The 
shaded areas correspond to the range 0.25 < |p|/T < 4 of 
momenta |p| of the decaying particle ipi with respect to the 
rest-frame of the medium. Here we assumed a thermal Bose- 
Einstein (BE) or Maxwell-Boltzmann (MB) distribution for 
b/b with vanishing chemical potential. In the low-temperature 
limit (NR), the vacuum value is approached. In the high- 
temperature limit (UR), the CP- violating parameter is en- 
hanced within the toy model. We also show the thermally 
averaged CP- violating parameter (ei) for the BE (red long- 
dashed line) and MB (blue dashed line) cases. 



As we have argued above, the conventional interpre- 
tation of the CP-violating parameter is only possible if 
the off-diagonal components of the full propagator peak 
at the same values of the momentum as the diagonal 
ones. In Fig. 4 we show the qualitative behavior of the 
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FIG. 4: Qualitative behavior of the components of the full 
and diagonal spectral functions and Q l p for Fj <C Mj and 



M 3 T 3 



|AMg| 



For illustration we choose M2/M1 = 2, 
gi\/M\ = 0.5 and I52I/M2 = 0.8. The diagonal components 
are down-scaled by a factor of 50. 



components of the diagonal and full spectral functions as 
obtained from Eq. (B50). As one can infer from the plot, 
the off-diagonal components G*J do peak at the same 
values of the momentum argument as the diagonal ones. 
Furthermore, the peaks of G" and Q™ are almost in- 
distinguishable. Therefore we can use the Boltzmann 
equations for the diagonal propagators Q to calculate the 
asymmetry. 

To conclude this section, let us estimate the range of 
applicability of (38). It was obtained by approximating 
the denominator of (27) by unity. Since Q\ = Gr and 
fr A = flu we obtain 



M 2 (Al- 



E 



[2\ 9j \ 2 \G 3 M\ 2 



(40) 



4Re (<?,<?* 



Evaluated on the mass-shell of the i'th quasiparticle, the 
retarded propagator, 

1 



Mf 



ll R 



(41) 



is of the order of the inverse splitting of the squared 
masses, I /AM 2 . Thus the contribution |C/^IT^| 2 is of 
order of (II/AM 2 ) 2 , whereas the contribution involving 
Re(^'n^) is of order of (11/ AM 2 ). 

In the hierarchical limit these ratios are much smaller 
than unity and the denominator of (27) can well be ap- 
proximated by 1. On the contrary, in the quasidegener- 
ate case the corrections can become large and (38) is no 
longer applicable. 



B. Resonant case 

If the spectrum of heavy (toy-) neutrinos is quaside- 
generate, \M\ — Mi\jM\ <C 1, the CP-violation parame- 
ter (2) predicts a resonant enhancement of the generated 
asymmetry, known as resonant leptogenesis [6]. Because 
of the enhancement, this scenario allows to circumvent 
the lower bound on the lightest right-handed neutrino 
mass Mi > 10 9 GeV typical for thermal leptogenesis [49]. 
Therefore, resonant leptogenesis is discussed as a possi- 
bility to evade constraints from the production of grav- 
itinos due to their impact on big bang nucleosynthesis, 
associated with the necessity of reheating temperatures 
well above Mi in the hierarchical case [50] . Furthermore, 
it has even been argued that the resonant enhancement 
could lower the scale of right-handed neutrinos to the 
TeV range, with possible implications for collider exper- 
iments [8, 51-53]. Therefore, it is important to check 
whether the conventional Boltzmann treatment which 
uses the canonical CP-violation parameter (2) agrees 
with the nonequilibrium field theory description in the 
quasidegenerate limit. 

Let us start by noting that the result (38) for the CP- 
violating parameter in the strongly hierarchical limit for- 
mally agrees with the canonical result (2) even in the 
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quasidcgenerate case. However, its derivation did involve 
approximations that break down in the resonant case. 
We will now revisit this derivation, starting from Eq. (27), 
without the above approximations. As explained, the de- 
nominator of (27) can significantly deviate from unity in 
the resonant case, which we take into account here. Using 
Eq. (40) we find 



e; = 



-2Im( 5j / 5i )lm(^n 



1 + 



|g.j| 2 

|9*I 2 



\g»n>*\*-2Tte(g J /g i )lte(g»W*) 



(42) 



Note again that we implicitly assume summation over j ^ 
i in the enumerator and the denominator, respectively. 

From Eqs. (40) and (32) we see that the denominator 
of (42) involves the self-energy IT^ , which is logarithmi- 
cally UV- divergent. It can be renormalized by includ- 
ing a mass counterterm. As shown in Appendix D, this 
amounts to the replacements Q — » Q ren , and 

tl R{ A)(X,p) -+ fl r R f A) (X,p) = n R{A) (x, P ) - 5Z P 2 + SM 2 . 

(43) 



For our purposes, it is convenient to use an on-shell renor- 
malization scheme, for which 



nr 



i/ 2 
(P 



nren.ij / 2 
h,vac(P 



M 2 ) 
M 2 ) 



(i = 1, 2) 



_n re ™' i Vn 2 - M 2 ) 

d p 2 LL h,vac \P ~ IV1 i ) 



(i = 1,2) 



(44) 



where W h ™ ac 



Ren^f™ ac denotes the dis- 



persive part of the renormalized self-energy in vacuum. 

Note that, in vacuum, the self-energy is time- 
independent and depends only on p 2 due to Lorentz in- 
variance (there is no medium which singles out a pre- 
ferred frame). Note also that the six independent coun- 
terterms described by the symmetric two-by-two matrices 
8Z and 5M 2 are fully determined by the above renor- 
malization conditions. For the renormalized self-energy 
in medium, we find (see Appendixes A and G) 



nr 



l (x, P ) 



+ 



l (X,p) 

,2 / | p 2| 



16tt 2 



hi- 



M? 



rr™(A, P ) = 



Rc(g t g|) 

16tt 2 

p 2 -m 2 

M 2 - M 2 



+ ■ 



+ 



Mf 



p 2 -m; 
m? 



- M 2 M? 



(45) 



In 



Mf 



(i ^ J) , 



where, in the medium rest-frame (for i, j = 1, 2), 
Re(^5* 



n™ ed ™(X,p) = 



x In 



16tt 2 |p| 

(2E- 



dE [f b (E) 



Pi 



(2E- 



Pl 



h(E)] 

(46) 



It is also straightforward to generalize the decomposi- 
tion (32) to the renormalized case. While the imaginary 
part L p (X,p) is not affected, we define (for i, j = 1,2) 



H™(A» ee ^{g i9 * + g ; gj )L™"(X,p) . (47) 

Since we will only use the renormalized quantities from 
now on, we omit the superscript 'ren' for brevity. Then, 
using also (36), we find 



16tt 



Im 



9*9j 



(48) 



3 3 \9i\ 2 \9j\ 2 



[Mjr i \0%L»\* + MS 3 M)] 



-1J 

J R 



T 13 



\L P . Note 



where MjTj = \g 3 ;| 2 /(167r) and L 
that , which is given by (41), as well as its real and 
imaginary components, Q" and Q 33 , also contain the 
renormalized self-energy. 

The propagators and the components of the self-energy 
in (48) are functions of the center coordinate X and 
the four-momentum p. In general, the components of 
p are not related by the on-shell condition, so that 
(48) describes CP-violating effects in the decays of on- 
and off-shell toy-Majoranas. However, if \M 2 — M 2 \ > 
MiTi,MjTj then the largest contribution to comes 
from on-shell momenta. The in-medium dispersion re- 
lation E™ ed is determined by the equation 4 



p 2 -M 2 ~ ReU%(X,p)\ p0=Erd = . 



(49) 



An effective, momentum-dependent mass can be defined 
by the relation E™ ed = [{Mf led ) 2 + |p| 2 ]i, which implies 



M, 



med 



(M 2 + ReU%(X,p)\ p „ =Erd y . 



(50) 



By means of Eq. (50) we take the mass-shift of the in- 
medium masses, 8 Mi = MJ 11 ^ — Mi, into account. In 



4 The contributions of off-shell momenta can be partially taken 
into account by evaluating the CP-violating parameter at a com- 
plex value of p determined by p 2 — M? — 11^ (X, p) = [24]. 
However, note that this approach (just as the original result (48), 
which is valid also off-shell) requires the use of an off-shell gen- 
eralization of the Boltzmann equations. 



9 



a thermal medium, they correspond to thermal masses. 
Even though the mass shifts 8Mi are small compared to 
Mi, they can be important for the mass difference in the 
quasidegenerate case. 

Finally we evaluate Eq. (48) for p° = E^ ned , using 
Eqs. (41) and (36). We obtain the following result for 
the CP-violating parameter: 



16tt 



Im 



9i9j 
9t9j 



(51) 



A, 



-"CP [Ay 



'CP 



[A 2 - 

L ij 



■ L 



p > 



encodes the medium cor- 



where L p = L p (X,p)\ p0=1 
rection discussed already in the hierarchical case, ob- 
tained by replacing Mi — > M™ ed in (35). The denom- 
inator consists of two contributions, weighted by the pa- 
rameters 



C CP 



[Re( fflff *)/(| 5i || 5 ,|)] : 



cos 2 (5 C p) 



(52) 



s 2 C p = [lra( gi g*) / (\ gi \\g j\)} 2 = sin 2 (6 C p) , 



related to the CP-violating phase Sep = arg^g^)- Fur- 
thermore, the degeneracy parameter, 



(M™ erf ) 2 - (M™ ed ) 2 + 5\ 



(53) 



is given by the difference of the in-medium ("thermal") 
masses of the heavy toy-neutrinos, plus logarithmic cor- 
rections (here we display the superscript 'ren' again for 
clarity): 



5'; 



2M ] T J L r r^(X lP )\ p0=Erd , 
2M J T J ( L r^{X,p)\ p0=ET , d 

-L 7 h en ' n {X,p)\ p0=E „i, d 



(54) 



As we have already mentioned, the resonance effects 
modify not only the CP-violating parameters, but also 
the total decay widths of the heavy particles; see Eq. (28). 
Using the functions and definitions introduced above we 
can write the in-medium decay widths in the form 



= r,x 



(55) 



°CP 



[A 2 - 



(^■r,) 2 ! 2 ] 



[Ai 



2c 2 7P ^f + [ S 2 7P (M,r j P p ) + cl P 5 2 J{M 3 T 3 L p )] 



In the limit \Mf - Mf\ > MjTj, the CP- violating pa- 
rameter (51) converges toward the result (38) for the hi- 
erarchical case. Our result for the resonant case, Eq. (51), 
describes the leading corrections when the above ratio be- 
comes sizable. However, one should keep in mind that, 
in the "maximal resonant" case, where the decay width 
is comparable to the mass difference, the Boltzmann pic- 
ture breaks down. The following simple argument sup- 
ports this statement: To a good approximation the spec- 
tral functions G" of the heavy fields have a Breit-Wigner 



shape. To evaluate the gain and loss terms (5) we have 
to integrate over the frequency k° of the heavy particles' 
propagator. If the distance between the peaks of the 
spectral function is considerably larger than the decay 
widths, the integration reduces to two independent inte- 
grations in the vicinities of the corresponding mass shells. 
Thus, we can identify two independent quasiparticle ex- 
citations with the corresponding distribution functions 
fip 1 and f^ 2 and the CP-violating parameters e\ and C2 
which contribute to the generation of the asymmetry On 
the other hand, if the decay width is comparable to the 
difference of the masses, then the peaks of the spectral 
functions overlap and neither the distribution functions, 
nor the CP-violating parameters, are well-defined. In 




FIG. 5: Qualitative behavior of the diagonal and off-diagonal 
components of full and diagonal spectral functions G p ] and 
g™ in the maximal resonant limit M] > \AMfj\ < MjTj. 
For illustration we choose M%/M\ = 1.02, |ffi|/Mi = 0.5 and 
1 02 1 /Ma = 0.8. The diagonal components are not scaled. 

Fig. 5 we present the qualitative behavior of the diago- 
nal and off-diagonal components of the full and diagonal 
spectral functions G l J and Q". As one can infer from 
the plot, the off-diagonal components of the full spectral 
function no longer have two pronounced peaks. Thus, 
it is not possible to define two CP-violating parameters. 
The shape of the diagonal components deviates from the 
Breit-Wigner one, and the positions of the peaks of G l p 
are shifted as compared to those for the diagonal spec- 
tral function Q l p . In other words, using the Boltzmann 

equation for the diagonal propagators Q to calculate the 
generated asymmetry is no longer a good approximation 
because of the off-shell effects. Moreover, since the micro- 
scopic time scales (t m i c ~ AM -1 ) and the macroscopic 
time scales (t mac ~ T _1 ) can be of the same order of 
magnitude in the "maximal resonant" regime, the mem- 
ory effects can play an important role. Thus we conclude 
that in this case, one should use the quantum kinetic 
equations discussed in Appendix B, which take off-shell 
effects into account, or even the full nonlocal quantum 
evolution equations discussed in Appendix E, which also 
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capture memory and correlation effects. 

With these considerations in mind, let us now ana- 
lyze the Kadanoff-Baym/Schwinger-Keldysh result (51) 
for the CP-violating parameter in the resonant (but not 
maximal resonant) regime, for some particular cases of 
interest. For definitcness we concentrate here on ei, i.e. 
i = l,j = 2 assuming that Mi < M 2 . In terms of the 
ratio 



R. 



Mf 



Mf 



(56) 



the condition of validity for the Boltzmann-treatment can 
be expressed as R 3> 1. The hierarchical limit corre- 
sponds to R — > oo. In the following we study the leading 
corrections in 1/R. 

In order to separate medium and resonance correc- 
tions, we first consider the vacuum limit /, 



6(6) 



0. 



The renormalization prescription Eqs. (44) ensures that 

M rned _^ M vac = M ^ ; & M[ &nd ^ are thc on _ shcll 

masses in vacuum. Furthermore, adopting the scheme 



(44) also implies that 5i 



S vac 



and 



\9jV_ 
16tt 2 




Mf 



Mf 



Mf 



(57) 

Therefore, in the vacuum limit the degeneracy param- 
eter is given by 



A*j -> &"ij C = Mf — Mf + {5[ 3 f 



(58) 



For \gj\/Mj <C 1, the loop correction can be safely 
neglected. Thus, in the vacuum limit the CP-violating 
parameter is approximately given by 



\9j\_ 
16tt 



-Im 



Hi. 



9i9j 



Mf - Mf 



(M 2 -M 2 ) 2 +S 2 cp M 2 T 2 



sm(25 C p) 



R 



R 2 + sin 2 (S C p) 



(59) 



There are two sources of medium corrections in the res- 
onant case. The first is the contribution from the spec- 
tral (imaginary) part of the self-energy loop, given by 
L p in Eq. (51). This contribution is known already from 
the hierarchical case; see Eq. (38) as well as Fig. 3. The 
second contribution stemming from the real part of the 
self-energy loop, given by L™ ed , enters via the thermal 
masses and also via the "logarithmic" corrections. 

For illustration, we insert a Bose-Einstein distribution 
for fb and /j (with \i = 0) in the expression for n™ ed in 
Eq. (46). In the limit T < M u we obtain 



(M™ ed ) 2 ss Mf 




l\9i\ 2 

6 TvT 3 " 



l I* I 2 

6 "MF 



(60) 



(M™ ed ) 2 - (M med ) 2 



Mf 



(M™ ed )- 2 ) T 2 . 



(61) 



The term in square brackets is of second order in the 
mass-squared difference and can be neglected. Assuming 
in addition \gj\ 2 <C Mf , the degeneracy parameter is 
approximately given by 



An ^Mf-M 2 + — 




Mf 



(62) 



Compared to the conventional result (2) the width in the 
denominator is effectively changed according to Tj — > 
\scp\Fj. Since \scp\ < 1, this means that for fixed val- 
ues of R and 5c p the result above is always larger com- 
pared to the conventional one. This property may be 
attributed to the fact that the Kadanoff-Baym formal- 
ism takes a resummation of resonant contributions into 
account. Note that, if taken at face value, the expression 
(59) for e\ ac formally has a peak at R„ ~ |scp|> with 
maximum |e*| = \ccp\ and width ~ jscplTj. The width g i rj sin(2Scp) < 
reduces to zero in the limit of vanishing CP-violation, 
5c p — > 0, although the peak value remains finite. On the 
contrary, for the conventional result (2) one has Pt* <~ 1, 
\e*\ = \scpccp\ and a width ~ Tj. However, we stress 
that the result should only be trusted in the regime where 
R » 1, and can be modified significantly by medium ef- 
fects. 



In the quasi degenerate limit it is plausible to assume that 
also the couplings are quasi degenerate, |gi| 2 ~ |ff2| 2 - In 
this case, the thermal corrections stemming from the dis- 
persive part of the self-energy cancel approximately. In 
general, the thermal corrections could also be larger than 
the mass-splitting of the vacuum masses. If the vacuum 
masses are quasi degenerate, this would destroy the reso- 
nance condition. However, also the opposite case is possi- 
ble, namely Ay could become tiny due to a cancellation 
of the vacuum masses and the thermal corrections for 
a certain temperature. Then the resonant enhancement 
would occur only close to this particular temperature. 
We do not pursue these possibilities further here. As- 
suming, for simplicity, that the thermal contributions in 
Eqs. (61) and (62) are subdominant, we can summarize 
the results for thc medium- and resonance corrections to 
the CP-violating parameter as follows: 



R 



R 2 + l 
R ■ L n 



R 2 
RL 



R 2 + . 



,2 T2 
CP L p 



conventional 



mod. corr. 



med.+res. corr. 



(63) 



The first two expressions follow directly from Eq. (2) and 
Eq. (38), respectively. The third one approximates the 
resonant result, Eq. (51), for R > 5. 
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In Fig. 6 we show the ratio of the expressions for 
the CP-violating parameter, Eq. (51), which includes 
medium and resonance corrections, and Eq. (38), which 
includes only medium corrections, for several values of 
the degeneracy parameter R. Large values of R corre- 
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FIG. 6: Corrections to the effective self-energy CP-violating 
parameter ei in the resonant case obtained from the Kada- 
noff-Baym formalism (for \gi\ = |ga| = O.OlMi, Sep = 7r/8 
and |p| = T; the mass M2 is determined by the value of R). 

spond to the hierarchical case and there is no difference 
between the two approximations. On the other hand, 
for small values of R the corrections are substantial. For 
instance for R = 10 the expression without resonance 
corrections underestimates the CP-violating parameter 
by a factor of two at high temperatures. 

Since typically most of the asymmetry is generated at 
T ~ Mi, it is instructive to look at the CP- violating 
parameter as a function of the degeneracy parameter R 
at T= Mi. 

In Fig. 7 we present the P-dependence of the con- 
ventional vacuum approximation for the CP-violating 
parameter, Eq. (2), the hierarchical approximation in 
medium, Eq. (38), and the resonant expression, Eq. (51), 
respectively. Very large values of R correspond to M 2 ^> 
Mi. In this case the resonance effects are suppressed 
and all three expressions go to zero oc P _1 . For smaller 
values of R we observe a significant deviation of the CP- 
violating parameter from its vacuum value, which is due 
to medium effects. Finally for even smaller values of R 
the resonant effects become important and we observe a 
deviation of the CP-violating parameter from its value 
calculated in the hierarchical approximation. 

The effective decay widths rj™ e ' 1 are also enhanced by 
the medium and resonance effects. The enhancement in- 
creases with the temperature; see Fig. 8. Furthermore, it 
strongly depends on the values of the degeneracy param- 
eter R, as can be seen in Fig. 9. Even for reasonable val- 
ues of the degeneracy parameter the effective decay width 
can be twice as large as in vacuum. This leads to a faster 



with med.+res. corr. 
with med. corr. 
conventional 




10 W z 

R = \M% - Ml\/M 2 r 2 



FIG. 7: Dependence of the vacuum (dotted line), hierarchi- 
cal (dashed line) and resonant (solid line) approximations for 
the CP-violating parameter on the degeneracy parameter R 
calculated at T — Mi (for the same parameter values as in 
Fig. 6). The dot-dashed line shows the approximate expres- 
sion for the resonant case from the last line of Eq. (63). The 
Boltzmann-approximation requires R S> 1 and therefore the 
results for the CP-violating parameter are not applicable in 
the gray shaded region. We show them only for comparison 
with the conventional result. 
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FIG. 8: Ratio of the in-medium decay width of the lightest 
toy-Majorana to its tree-level value as a function of the di- 
mensionless inverse temperature calculated for various values 
of the degeneracy parameter R. 



decay of the toy-Majoranas. Also the inverse decay pro- 
cesses are more efficient in washing out the asymmetry 
In other words, the increase in the generated asymmetry 
due to the enhancement of the CP- violating parameter 
can be partially compensated by the increase of the in- 
medium decay widths. Let us also note that the observed 
enhancement of the total in-medium decay widths could 
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FIG. 9: Ratio of the in-medium decay width of the lightest 
toy-Majorana to its tree-level value as a function of the de- 
generacy parameter R calculated at T — Mi . 

be very important in the scenarios with very small values 
of the degeneracy parameter [54, 55]. 

IV. NUMERICAL RESULTS 

For the moment we consider numerical solutions only 
for the strongly hierarchical case, as the resonant case is 



more involved and one cannot obtain consistent Boltz- 
mann equations for the maximal resonant case. To ob- 
tain the Boltzmann equations for f b and f b we integrate 
Eqs. (3) and the corresponding equation for 6, together 
with Eq. (5), over po- The Boltzmann equations for f^ )i 
are obtained from Eq. (B52). As one can infer from Fig. 4, 
in the hierarchical case the off-diagonal components of 
the full propagators G are subdominant and the diag- 
onal components of G are almost identical to those of 
Q. Therefore, we can neglect the off-diagonal compo- 
nents in the kinetic equations (B52) for the full propa- 
gators and approximate them by the kinetic equa- 
tions (B53) for the corresponding diagonal propagators 
Q^. The Boltzmann equations are then obtained after 
using the Kadanoff-Baym ansatz (8) and the quasiparti- 
cle approximation (9). 

We solve the coupled system of Boltzmann equa- 
tions in the spatially homogeneous and isotropic case 
in (spatially flat and radiation dominated) Friedman- 
Robertson-Walker space-time. They take the form 

L[/](|p|)^2? Q /(|p|) 

= P °(|-' P '^) /(IPl) ' (64) 

where H = a /a is the Hubble parameter. As usual, the 
integrations over the time components of each of the in- 
variant four-volume elements in the collision terms can 
be performed trivially after the quasiparticle approxima- 
tions for the spectral functions have been inserted. 



The resulting system of Boltzmann equations takes the same form as the one presented in [2] for the vertex 
contributions. As in the vertex case the structure differs from the usual one obtained in the conventional bottom-up 
approach. In particular, we do not need to include the RIS part of the collision terms for the processes bb <-> bb 
because our collision terms for the processes bb O ipi and bb O ipi do not suffer from the generation of an asymmetry 
in equilibrium. The form of these equations is also necessary to guarantee cancellation of the gain- and loss-term 
contributions in equilibrium when the quantum statistical terms are present. This structure can be translated directly 
from the toy-model to established scenarios of leptogenesis and baryogenesis by analogy. Therefore, we consider it as 
important, also for phenomenological studies, and repeat it here: 

L[A](|P|) = C , b 6^6[/b](|p|)+C' 6 - b ^[/ b ,/5](|p|) + C 6b ^ 1 [/ b ,/^](|p|) , (65a) 
L{M\P\) = C 55o55 [/ 5 ](|p|)+C 5bo5b [/ 5 ,/ 6 ](|p|) + C 55 ^ 1 [/ 5 ,/ v , 1 ](|p|) , (65b) 
L [JW ] (I p| ) = « bb [U ± , f b ] (|p|) + ^ bb [U t , f b ] (|p|) . (65c) 

Here, Cj<->/ denotes the collision term for a process i — > f. The collision terms for the 2 — 2 scattering processes in 
(65a) are given by 

[A] (IpI) - Udn 3 P2 dn 3 P3 dn 3 P4 (2TT) i s(p + P2 - P3 - Pi )\\ 2 

x{[l + A(|p|)][l + /6(|p 2 |)]/6(|P3|)/6(|P4|) - A(|p|)/6(|pa|)[l + A(|Ps|)][l + A(|P4|)]} , (66a) 

c b - b ^ b - b [hJ- b ](\p\) = ydii 3 P2 dn 3 P3 dn 3 P4 (2TT) i s(p + p 2 ~ P3 - P4 )A 2 

X{[1 + A(|p|)][l + A(| P2 |)]A(|P3|)/ 6 (|P4|) - /6(|P|)/6(|P 2 |)[1 + /6(|P3|)][1 + A(|P4|)]} ■ (66b) 

Replacing f b with in Eqs. (66) one obtains the analogous terms in the equation for b. The collision terms in Eq. (65c) 
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are obtained by inserting the diagonal components of the self-energy Eq. (F5) into Eq. (B53): 

~ \ JdH 3 p2 dU 3 p3 (2TT) 4 5{p-p2-p 3 )±\ gi \ 2 

X [{[1 + jW(|p|)]/ 6 (|P2|)/6(|p3|) - /^(| P |)[1 + A(|p 2 |)][l + /6(|P3|)]} 

+ {[1 + /^(IPI)]/5(IP2|)/6(|P3|) - /W(|p|)[l + h(\P2\W + h(\Ps\)}}] ■ (67) 

In the framework of the toy model the CP-violating parameter for the self-energy loop contributions given in 
Eq. (38), in the strongly hierarchical limit AI2 ^ Mi, differs by just a (symmetrization) factor 1/2 from the vertex 
contributions. It appears explicitly in the collision terms for the (inverse) decay of ip\ into bb or bb: 

C , 6fc «.0j/ 6 ,/ v J(|p|) = i/dn3 2 dn3 3 (27r) 4 ( 5(p 2 -p-p 3 )| 3l | 2 [l + ei(|p3|)] 

x{[l + /„(|p|)][l + f b (\p 2 \)}UA\P3\) - A(|P|)/6(|P2|)[1 + UM)}} , (68a) 
^^J/^^ 1 ](|p|) = H dn p 2 dn p 3 ( 2 ^fe"P-P3)l3i| 2 [l-ei(|p 3 |^ 

x{[l + / 5 (|p|)][l + / 5 (| P2 |)]/^(|P3|) - /b(|p|)/6(|P2|)[l + /^(|P3|)]} • (68b) 



The network of Boltzmann equations (65) should be un- 
derstood in a generalized sense. The "amplitudes" which 
appear here differ from the usual perturbative matrix el- 
ements and do not share their symmetry properties. 

In order to study the effect of the quantum correc- 
tions, we can again compare the results obtained by 
integrating the network of Boltzmann equations with 
quantum-corrected ei(|p|) with the corresponding ones 
in the vacuum limit e\ ac . The computation is started 
at sufficiently high temperatures so that all species, in- 
cluding -01 with mass M\ = 10 10 GeV, have relativistic 
initial distributions. In addition, we assume that the in- 
teractions are in chemical equilibrium in the beginning, 
i.e. fj,^ 1 = 2/ib = ^P"b- We start with sufficiently neg- 
ative chemical potentials fib as to avoid Bose-Einstein 
condensation of the different species 5 . 

We choose the coupling A (and gi via ei and k) such 
that the 2 — 2 scattering rates are much larger than those 
of the decays and inverse decays. This assures that the 
light species are kept in kinetic equilibrium, as in the 
standard leptogenesis scenario. As shown in [2] there is 
no need to compute the collision integrals for 2 — 2 scat- 
tering explicitly in this case. This means that they can 
be described in terms of four parameters fib, 7& and fib, 
Xg which obey the relation Tj = Tf,. Hence, we studied 
the evolution of fb and fb in terms of only three parame- 
ters. In contrast, the full equation for -0i was discretized 



As was explained in [2] this necessity arises here, because we 
consider the system (65)-(68) as closed so that there are no in- 
teractions which can remove the produced over-densities of fe's 
and fe's from the system. Therefore 6 and b can in principle un- 
dergo Bosc— Einstein condensation, which we avoid by choosing 
IH, and fif, appropriately. Such interactions will be present in 
a phenomenological scenario. Whether the possibility of Bose- 
Einstein condensation exists in such scenarios will have to be 
answered by solving appropriate kinetic equations. 
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FIG. 10: The distribution function for the massive species fy 1 
can deviate significantly from equilibrium for smaller washout 
factors (here k, ~ 0.01). 



on a grid with 400 momentum modes and solved simul- 
taneously. For this purpose, homogeneity and isotropy 
was assumed, so that the angular integration can be per- 
formed as in [56] in a rather general case or in the ap- 
pendix of [2] for the present special case. As shown in 
Fig. 10, for small washout parameters, 

K = Ti/H{M X ) = | 5l | 2 m pi /(4.5 • 16^ V^?) , (69) 

the distribution function can deviate significantly from 
the equilibrium form (for which the curves would be 
straight lines). An equilibrium form would be a neces- 
sary assumption to obtain rate equations. The generated 
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FIG. 11: The ratio (ei) /ei ac - The curves flatten for small 
Mi /T because the initial conditions involve a finite chemical 
potential. 



FIG. 12: The final asymmetries and the relative quantum 
correction {r] — rf ac )/rf ac in dependence of the washout factor 
k. The cases a, b, c, d, e, f correspond to washout factors 
0.01, 0.1, 0.366, 1, 10, 100. 



"baryon" asymmetry is defined as 

,(«,/^=»S2, (70, 
s{Mx/T) 

where rtf, and are the number densities of species b and 
b and s is the standard cosmological entropy density [25] . 
We denote the analogous quantity, corresponding to the 
solution for e\ ac , by rf ac (Mi/T). Fig. 11 shows the result 
for the ratio (ei) /e\ ac for various values of the washout 
parameter. Comparing it to the thermal equilibrium re- 
sult in Fig. 3 one sees a flattening for small M\jT which 
is caused by the finite chemical potential of b and b in the 
initial conditions. One would obtain larger corrections if 
additional interactions for b and b would be introduced 
in order to start with smaller chemical potentials. 

The generated "baryon" asymmetry does not depend 
monotonously on the washout parameter k when the 
medium corrections are taken into account. This can 
be inferred from Fig. 12 where the dependence of the 
final asymmetries r/ = rj(Mi/T —> 00) and rf ac = 
r\ vac {M\jT ~ > 00) are presented. In the present case, 
where b and b are bosons, the quantum corrections always 
lead to an enhancement of the asymmetry compared to 
the results without the corrections and the asymmetry 77 
has a maximum for moderate washout factors k ~ 0.059. 
The maximum of the relative enhancement of about 26% 
is reached at k ~ 0.34. The enhancement of t\ due to 
the quantum corrections is suppressed at large washout 
factors, since the same processes which create and dimin- 
ish the asymmetry, are effective at late times, where the 
CP-violating parameter takes smaller values (compare 
Fig. 11). In the opposite case of small k, the particles de- 
cay late so that the washout is ineffective. However, the 
interval of integration in Eq. (35) is located at relatively 
large momenta since the mass increasingly dominates the 



relativistic energies as the momenta are red shifted to 
smaller values. This means that the integration is over 
an interval in which the distribution becomes smaller 
and smaller. Therefore, the relative quantum corrections 
go to zero for small k. 

This interpretation of the results was already given in 
[2] as they are the same for the vertex and self-energy con- 
tributions in the hierarchical limit. We can here draw the 
additional conclusion that the combined effect from both 
contributions is the same in this limit. This is important 
since they could in principle have opposite effects. This 
is not the case as the self-energy and vertex CP-violating 
parameters differ only by a positive prefactor, just as in 
the vacuum case. However, differences between the two 
can appear in the resonant regime. In this case the ex- 
pressions for the vertex and self-energy contributions to 
ei have a different momentum dependence. Additionally, 
in general one has to take into account a further Boltz- 
mann equation for the second heavy species ip2, which is 
currently being studied. We would like to stress again 
that the size and sign of the corrections depend on the 
quantum statistics of the particles in the vertex- and self- 
energy loops and will be different in a phenomenological 
scenario. 



V. CONCLUSIONS AND OUTLOOK 

In this paper, we have studied leptogenesis in a simple 
toy model consisting of one complex and two real scalar 
fields in a top-down approach, using the Schwinger- 
Keldysh/Kadanoff-Baym formalism as starting point. 
This treatment, based on nonequilibrium quantum field 
theory techniques, is motivated by the fact that it al- 
lows a unified description of two key ingredients of lepto- 



15 



genesis, namely deviation from thermal equilibrium and 
loop-induced CP -violation. 

We find that the structure of the kinetic equations au- 
tomatically ensures that no asymmetry is produced in 
thermal equilibrium. In other words there is no need for 
the real intermediate state subtraction, i.e. the formal- 
ism is free of the double-counting problem typical for the 
canonical approach. 

One of the key quantities in leptogenesis is the CP- 
violating parameter. Earlier studies have shown that 
there are two sources of CP-violation, the vertex and the 
self-energy contribution. In this work, we have concen- 
trated on the latter. We have found that for scalar fields 
medium effects increase the self-energy contribution to 
the CP-violating parameter. 

Contrary to the results obtained earlier in the frame- 
work of thermal field theory by replacing the zero tem- 
perature propagators with finite temperature propaga- 
tors in the matrix elements of the Boltzmann equation, 
the medium corrections depend only linearly on the par- 
ticle number densities. 

Although the formal description of the self-energy and 
vertex contributions to the CP-violating parameter is 
technically quite different, the results for both are very 
similar qualitatively and, in the hierarchical case, even 
quantitatively. In this work, we additionally studied the 
quasidegenerate case, for which the self-energy contribu- 
tion is essential. 

We have shown that the canonical expression for the 
self-energy CP-violating parameter is only applicable in 
the hierarchical case, even though it docs not diverge in 
the limit of equal masses. In the resonant regime the 
interactions modify the mass spectrum of the quasiparti- 
cle excitations. Furthermore, using the Kadanoff-Baym 
formalism, it is possible to take a resummation of res- 
onant contributions into account. Both effects lead to 
changes in the expression for the CP-violating parame- 
ter. For moderate values of the degeneracy parameter R 
the resonance corrections can enhance the CP-violating 
parameter by a factor of two. Therefore, it is impor- 
tant to take these corrections into account in numerical 
simulations. 

Another important effect is the resonant and medium 
enhancement of the total decay widths. It leads to a 
faster decay of the heavy particles and more efficient 
washout of the generated asymmetry. Therefore, the in- 
crease in the generated asymmetry due to the enhance- 
ment of the CP- violating parameter can be partially com- 
pensated by the increase of the in-medium decay widths. 

In the "maximal resonant" regime the Boltzmann pic- 
ture is no longer applicable. This can be attributed to the 
fact, that in this regime the peaks of the spectral func- 
tions of the heavy (toy-)neutrino fields overlap and it is 
no longer possible to unambiguously define quasiparticles 
and one-particle distribution functions. Furthermore, the 
off-diagonal components of the correlation functions no 
longer have two pronounced peaks and therefore it is not 
possible to describe the mixing effects in terms of the 



corresponding on-shell CP-violating parameters. Conse- 
quently, in this regime the calculation of the generated 
asymmetry requires us to use at least two-by-two matrix 
equations for the diagonal and off-diagonal components 
of the propagators of the heavy fields with (in general) 
off-shell momenta. These can be obtained from the quan- 
tum kinetic equations by dropping the Poisson brackets 
on the right-hand side. Furthermore, since the micro- 
scopic time scales (i m , c ~ AM' 1 ) and the macroscopic 
time scales (t mac <~ T _1 ) can be of the same order of 
magnitude in the maximal resonant regime, memory ef- 
fects can play an important role. Taking both off-shell 
and memory effects into account consistently requires the 
use of the full system of Kadanoff-Baym equations. 

The formalism developed in this paper also provides 
a powerful tool for analyzing quantum noncquilibrium 
effects induced by the expansion of the early Universe. 
In particular, there is a small additional "spontaneous" 
contribution to the CP-violation in the system similar 
to that encountered in clectroweak baryogenesis [57-60] 
This effect will be investigated in a forthcoming paper 
[61]. 
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Appendix A: CP-violating parameter 

In this appendix, we review the calculation of the self- 
energy contribution to the CP-violating parameter in 
vacuum, e™ c = (r^_, 66 - I^.^) / (T^ bb + I^.^), 
in the conventional in-out formalism. 

Since the toy-Majoranas are unstable, they cannot ap- 
pear as in- or out-states of S-matrix elements. Instead, 
their properties are defined by S-matrix elements for scat- 
tering of stable particles mediated by the unstable neu- 
trino [62]. Resumming the propagator of the interme- 
diate heavy state, we can separate two-body scattering 
processes in resonance contributions and the rest. The 
CP- violating part of the resonance contribution can then 
be interpreted as a characteristic of the on-shell inter- 
mediate toy-Majorana [5, 63]. The amplitude of the s- 
channel two-body scattering process bb — > bb (see Fig. 13) 
can conveniently be expressed as 

AW 6 6 = £rfG«( s )rf , (Ai) 

where Tf and represent the vertices tpibb and ipjbb 
that include the wave functions of the initial and final 
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FIG. 13: Resummed diagrams contributing to the resonant 
part of the 2— »2 scattering amplitude of the process bb — > bb. 



states, and G lJ are the full propagators obtained by re- 
summing an infinite series of toy-Majorana self-energy 
graphs [5]. 

The resummation can be performed using the 
Schwingcr-Dyson equation in vacuum: 



[G- 1 ]«(p 2 ) = [ P 2 -Mf] ^-ir> 2 ) 

At one-loop level the self-energy IP-? reads 

n'V) = - g ' g L + f gj j V)> 



where 



32tt 2 



\p 2 



B(p z ) = A + 2 - In i±V + «T0 (P ) . 



(A2) 



(A3) 



(A4) 



7 + 4-7T contains the divergent contribu- 
tion. To renormalize the mass and the self-energy we 
introduce the wave-function and mass counterterms to 
the Lagrangian: 



6C 



1 



.V/.,,n> , r,i) i ,i- J - -8M?M 



(A5) 



where SZij and SM 2 j are symmetric two-by-two matrices. 
This implies that the renormalized self-energy is given by 



Kin(p 2 ) = n*V) - p 2 sz %3 + 8Mi 3 



(A6) 



In the on-shell renormalization scheme the dispersive 
parts of the components of the renormalized self-energy 
must satisfy the conditions 



Inverting (A2) we obtain for the components of the renor- 
malized resummed propagator: 

G u (p 2 ) = +[G-^V)/det[G- V)] , (A9a) 
G l \p 2 ) = -[G- 1 ] ii (p*)/det[G- l (p')] . (A9b) 

Because of the presence of absorptive terms in (A8) the 
determinant of the inverse propagator in (A9) has two 
poles in the complex plane at 



Mf - iMiTi , 



(A10) 



where T; = \gi\ 2 /WirAli is the tree-level decay width of 
ipi. Expanding (A9) around the poles and substituting 
the leading expansion terms to (Al) we find [6] 



M 



bb^bb 



"£vf(s)—V*(s), (All) 



where 

= r 



A(B) _ [G^fjs) 
[G-i]»(s) 

MB) ^ ir'» 



s-M?-n«(s) 



(A12) 



The modulo squared of the bb — > bb scattering amplitude 
is then given by 



\M 



bb^bbl 



E 



I^OOIWOOI 



MiTi 



cross terms. 



(s - Mi) 2 + (MiTi) 2 
(A13) 



The Breit-Wigner propagators in the diagonal terms of 
(A13) strongly peak on the mass shell of the quasipar- 
ticles, i.e. at s = M 2 and rapidly decrease off the mass 
shell. In the limit of vanishing decay widths 



{s - Mi) 2 + (M,r,) 2 



nS(s~M 2 ). 



(A14) 



IP* (v = M 



f2\ 



— IT" (v 2 = M 2 ) = 



(A7a) 



IP r i n (p 2 = M 2 )=H? en (p 2 =M 2 )=0 (i^j). (A7b) 

Using (A6) and (A7) and the explicit form of the bare 
self-energy we calculate SZij and <5M 2 . Substituting 
them into (A6) we find 



ren ^ 



111 



M 2 



P 



M 



M? 



T0(p 2 ) 



(A8a) 



n y = 

ren 



R*(ffifl?) 



16tt 2 



P 2 -M 2 J^| 
M 2 - M 2 M 2 

J J 



P~ M I ]n \P 2 \ 
M? - M? M? 



in6(p 2 



(A8b) 



Furthermore, if \M 2 - M 2 \ > ,. , then the two 
Breit-Wigner propagators do not overlap and we can ne- 
glect the cross terms in (A13). In other words, in this ap- 
proximation the resonant (real intermediate state) part 
of the amplitude is given by 



\M 



bb^bbl 



E M|pP T , (s _ M . ) , (A15) 



MiTi 



Equation (A15) suggests, that should also be eval- 

uated at s = Mf. As follows from (A12) to leading order 
in the couplings its modulo squared can be represented 
in the form 



r?(i-£i 



(A16) 
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Using the explicit form of the renormalized self-energy 
(A8) we obtain for the CP-violating parameters in vac- 
uum: 



\g£ (g i9 *\ Mf-Mj 

167T m V^J(M|-Mf)2 + M|r2- 



(A17) 



Because of the presence of the MjTj term in the denomi- 
nator of (A17), it does not diverge if Mj — » Mj. However, 
the condition |M?-M?| > M^j, MjFj is not satisfied in 
this case and therefore the approximations we have made 
to derive (A17) are not valid in this limit. Let us also 
note that due to approximation (A14) the CP-violating 
parameter (A17) characterizes on-shell heavy particles. 

Integrals of the left- and right-hand sides of (A14) over 
s are equal only if we integrate in the range from — oo 
to +00. But in the s— channel the momentum trans- 
fer squared is always positive. Moreover, the functions 
V^ 3 ^ also depend on s. Thus, in the transition from 
(A13) to (A15), i.e. in replacing the Breit-Wigner prop- 
agator by a Dirac <5-function and evaluating V^ 3 " 1 at 
s = Mf we have made an approximation which, strictly 
speaking, is only valid in the limit I\ — > 0. We could 
perform the integration more carefully and take the fi- 
nite widths into account using Cauchy's integral theo- 
rem and evaluating V^ 3 ^ at the poles Sj = Mf — iMiTi. 
This leads to a slightly different expression for the CP- 
violating parameter [24]: 



\9i\ 



■Im 



Mf 



16* ^gJiMf-Mfy + iMjTj-M^- 

(A18) 

Equation (A18) can be considered as a better estimate 
of the CP-violating effects in the system. However, since 
the on-shell approximation (A14) no longer applies, one 
can not interpret (A18) as an expression for the CP- 
violating parameter of an on-shell state. In other words, 
strictly speaking, (A18) would require us to use an off- 
shell generalization of the Boltzmann equation. 



Let us also note that in the case \Mf 



Mi 



> 



MiTi, MjTj, which we have considered here, the differ- 
ence between the two expressions for the CP-violating 
parameter can be neglected. 



Appendix B: Mixing real scalar fields 

In this appendix, we derive the Kadanoff-Baym, quan- 
tum kinetic and Boltzmann equations for a system of two 
(or, in general, n) mixing real scalar fields. 



1. Schwinger— Dyson equation 

The generating functional of such a system reads 
Z = f SfoSfa exp[t(S + Jiil>i + ^iKij^j)] . (Bl) 



The bilinear external source is a n-by-n matrix with the 
property Kij(x,y) — Kji(y,x), where n — 2 in the toy- 
model. Note that the field and the external sources are 
defined on the positive and negative branches of a closed 
real-time contour; see Fig. 14. Throughout this work, we 
use the compact notation of Rcf. [29] , which avoids the 
doubling of the degrees of freedom. In particular, note 
that the indices i,j € {l,...,n} refer to the two real 
scalar fields tpi in the toy model. 



tinit 



FIG. 14: Closed real-time path C. 



The functional derivatives of the generating functional 
for connected Green's functions, W = —ilnZ, with re- 
spect to the external sources read 



dW 

dJi(x) 
dW 
dKij(x,y) 



= , (B2a) 
= $[G ji (y,x) + 9 i (x)Hf j (v)], (B2b) 



where = (tpi(x)) is the expectation value of the 

field operator. The propagator Gji is a n-by-n matrix; 
its off-diagonal components describe the mixing of the 
two fields. We emphasize again that the indices refer 
to the field content, and have nothing to do with the 
branches of the time-contour. 

Performing a Legendre transform of the generating 
functional for connected Green's functions, we obtain the 
effective action 

r = W - J,* 4 - \tv\KijGji] - l^iK^j . (B3) 

Its functional derivatives with respect to the expectation 
value and the propagator reproduce the external sources: 



ST 



ST 

5Gij(x,y) 



-Ji(x) - J @ 4 z K l3 {x, z)^j(z) , (B4a) 
-faifax), (B4b) 



where £^ 4 z = y^gd 4 z and g is the determinant of the 
space-time metric. 

Shifting the fields by their expectation values Vi - ► 
tpi + we can rewrite the effective action in the form 

r = - iln / ^ViW>2 exp[i(S* + JiVi + ^ipiKijipj)] 
+ S cl m-\tT[K ij Gj i ]. (B5) 

Next, we tentatively write the effective action in the form 

T ee S c ,m + § lndct [G- 1 ] + itr^/G^] + T 2 , (B6) 
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which defines the functional IV 

Differentiation of the third term on the right-hand side 
with respect to the field propagators yields the inverse 
free propagators. Note that the latter one is diagonal - 
in the absence of interactions the two fields do not mix: 



■ 1 (x,y)=i(n x + Mf)5 9 (x,y)S.. 



(B7) 



where S 9 (x, y) = (—g x )~^5(x, y)(—g y )~ 1L is a generalized 
Dirac <5-function. 

The second term on the right-hand side is defined by 
the path integral 



det 



G~ 



2ir 



exp tyj GyVj) ■ 



To calculate the functional derivative of In det [G 1 ] , we 
take into account that 



/ S' 4 zG^ l 1 k (u,z)G kn {z,v) = 6 9 (u 7 v)S mn . 
After some algebra and use of (B8), we obtain 
5 



SGij(x,y) 



In det [G- 1 ]=-Gj?{y,x). 



(B8) 



(B9) 



The functional derivative of (B6) with respect to G then 
reads 

= -\K ji (y,x). (BIO) 

The considered physical situation corresponds to vanish- 
ing sources. Introducing the self-energy, 



Uij(x,y) = 2i- 



ST 2 



5Gji(y,x) 

we can then rewrite (BIO) in the form 

Gr j 1 (x,y)=^ 1 (x,y)-U ij (x,y). 



(Bll) 



(B12) 



As can be inferred from (B12), off-diagonal components 
of Gij are induced by off-diagonal components of II jj . In 
the considered model they arise because both real scalars 
couple to the "baryons" . 



2. Kadanoff— Baym equations 

Convolving the Schwinger-Dyson equations (B12) with 
G from the right and using (B8), we obtain 

ip x + Mf]G ij (x,y) = 6 9 (x,y)6 ij 

+ f@ 4 zTL lk (x,z)G k3 (z iy ). (B13) 

The statistical propagators and spectral functions are 
defined by 

G%{x,y) = \{[il H {x)^ j {y)]+), (B14a) 
G i p \x,y)=i{[^ i (x),^ j (y)}.). (B14b) 



From the definitions (B14) it follows that 

G%(x,y) = G%(y,x), &> (x,y) = -&> {y,x) . (B15) 

The time-ordered Schwinger-Keldysh propagator (which 
is the analogue of the Feynman propagator on the closed 
time path) is a linear combination of these functions: 

G l0 {x lV ) = G%{x,y) - IsigncOr - y°)G p j (x,y). (B16) 

Upon use of the signum- and (5-function differentiation 
rules, the action of the O x operator on the second term in 
(B16) yields a product of g 00 <5(x°, y°) and V%G pij (x,y). 
Using the definition (B14) and the canonical commuta- 
tion relations [64] 

lim [V>i(z°,x),7rj(a: ,y)]_ =iS(x,y)S ij , (B17) 

y u — >x u 

where nj = .g 00 ^/ r g VoV'j ; we nn d f° r the derivative of 
the spectral function 



V x G i j(x,y) = 



(B18) 



Multiplication of (B18) by g 00 5(x°, y°) then gives the 
generalized Dirac ^-function S 9 (x, j/)<5 lJ : , which cancels 
the 5- function on the right-hand side of (B13). Further- 
more, we decompose the self-energy according to 

IMz, j/) = iT/fotf) - isign c (x° - y°)IP p j (x,y) . (B19) 

The resulting system of Kadanoff-Baym equations reads 
(see [40, 65] for more details): 



[n x + Mf]G%(x,y) = f$> 4 zir>?(x,z)G k p i(z,y) 
o 



J$> 4 zW p k (x,z)G k F j (z,y), (B20a) 



[n x + M^(x,y) = J^ 4 zUf(x,z)G k /(z,y). (B20b) 

In the limit of just one scalar field it reverts to the system 
derived in [19]. Numerical solutions of full Kadanoff 
Baym equations, for systems involving effectively a single 
degree of freedom in the scalar sector, have been studied 
e.g. in [13, 14, 66-72]. 

3. Quantum kinetics 

The system of Kadanoff-Baym equations can be 
rewritten in terms of the advanced and retarded prop- 
agators, G^ and G^: 

\n x + Mf]G% (p) {x, y ) = -!^ze{z°) 

x \n$ [p) {x,z)G k J{z,y) + T^{x,z)G k J (p) (z,y)]. (B21) 
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In order to close the system, Eqs. (B21) must be sup- 
plemented with the analogous equations for the retarded 
and advanced propagators: 

[U x + Mf]G% {A) {x,y)=5^x,y)^ 

-J^zH% A) (x,z)G% A) (z,y). (B22) 

From the definitions of the retarded and advanced 
propagators and relations (B15) one can infer that 

G i l(x,y) = e(x°-y )G i i(x,y) 

= -6{x° y°)Gf(y,x) = &2(y,x) . (B23) 

Therefore, after interchanging x and y and then i and j 
in (B21), we obtain 

[D y + Mf}G F{p) (x,y) = -J^z9(z°) 
x [(%(x,z)n k j w (z,y) + (% M {x,z)n k j(z,y)]. (B24) 

The same operation applied to (B22) yields 

[□„ + M]]G A{R)i . {x, y) = 6 9 (x, y)S t3 

- J 9^z G A[R)ik (x, z)U A{R)k . (z, y) . (B25) 

Following the usual procedure, we introduce the cen- 
ter and relative coordinates X and s [73]. The Wigner 
transform of the statistical propagator is defined by 

G%{X 1 p) = J—g x jd i S e^ s G%{X lS ). (B26) 

The definition of the Wigner transform of G Z J (X, s) dif- 
fers from (B26) by an overall factor of — i. From (B15) 
it then follows that 

G* F {X, p) = Gp (X, p) , G* p {X, p) = G T p {X, p) , (B27) 

where the hats denote the corresponding matrices and 
the superscript 'T' denotes transposition. As we see, the 
Wigner transforms are Hcrmitian matrices. The Wigner 
transforms of the retarded and advanced propagators are 
defined analogously to (B26). Using (B23) one can then 
show that 

G R (X,p) = G A (X,p). (B28) 

Just as in the case of a single real scalar field [19], from 
the definitions of the Wigner transforms, relation (B23) 
and the equality 9(s°) + 6(-s°) = 1 it follows that 

G R (X,p)-G A (X,p) = iG p (X,p). (B29) 

We could have derived Eqs. (B28) and (B29) using the 
spectral representation of the retarded and advanced 
propagators 

G R{A) (X, P ) = - r ^ p (X,u P ) _ 
y ' J 2n p — to ± le 



The spectral representation also implies that G R and G A 
can be represented as a linear combination of two Hcr- 
mitian matrices: 

G R(A) (X,p) - G h (X,p) ± § G p (X,p) . (B31) 

Let us now subtract (B24) from (B21) and Wigner trans- 
form the left- and the right-hand sides of the resulting 
equation. On the right-hand side we neglect the mem- 
ory effects, that is, we discard the function and re- 
place X xz , X zy by X xy . Furthermore, we keep only terms 
at most linear in the Wigner transform of the covariant 
derivative T> a ; see [19] for more details. Introducing 
the Poisson brackets, 

{A(X,p),B(X,p)} PB = —A(X,p)V a B(X,p) 

Op a 

-V a A(X,p)^-B(X,p), (B32) 

Op a 

the generalized Poisson brackets 6 

[A(X,p),B(X,p)} PB = -t[A(X,p),B(X,p)}^ (B33) 
+ \{A{X,p), B(X,p)} PB + \{A{X,p), B(X,p)} PB , 

and 

(X, P)=P 2 - M 2 - fl R{A) (X, p) (B34) 

we can rewrite the result in a compact form: 

[Cl R (X,p),G Hp) (X,p)] pB 

= G F{p) (X,p)n p (X,p) - fl F{p) (X,p)G p (X,p) 
-[G R (X,p),n F(p) (X,p)] pB . (B35) 

In the case of a single scalar field the commutator in 
(B33) vanishes and, using (B31), we can then show that 
the quantum kinetic equation (B35) for the spectral func- 
tion and statistical propagator revert to those derived in 
[19]. After some algebra we rewrite the quantum kinetic 
equation for the spectral function in the form: 

[n(x,p),G p (x,p)] PB 

= -[G h (X,p),fl p (X,p)} PB , (B36) 

where 

n = p 2 - M 2 - Itfc . (B37) 

The quantum kinetic equations must be supplemented by 
the corresponding constraint equations. To derive these, 
we add up the Wigner-transforms of (B24) and (B21). 



If the matrices A(X,p) and B(X,p) are hermitian and commute, 
then [A(X,p),B(X,p)] PB = {A(X, p), B(X, p)} PB . 
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To linear order in the gradients the resulting constraint 
equations read 

[Cl R (X,p),G F{p) (X,p)} e = [Il F(p) (X,p),G A (X,p)} @ 
+ i{fl F{p) (X,p),G A (X,p)} PB(B 
+ i{TL R (X,p),G F{p) (X,p)} Pm , (B38) 

where we have introduced [A, B]® = AB + W A^ and 
{A,B} PB(B = {A,B} PB + {B\ti} PB to shorten the 
notation. In the case of a single scalar field (B38) revert 
to those derived in [19]. After some algebra we can again 
simplify the constraint equation for the spectral function: 

[Cl(X,p),G p (X,p)U = [fl p (X,p),G A (X,p)] (B 

+ i{fi p (x, P ),G h (x, P )} Pm 

+ i{fl h (X,p),G p (X,p)} Pm . (B39) 

To close the system of the constraint equations (B39) we 
have to derive the constraint equation for the function 
Gh or, alternatively, for the retarded propagator: 

[Cl R (X,p),G R (X,p)} + = -2I 

+ ±{fl R (X,p),G R (X,p)} PB+ , (B40) 

where the + operation is defined as the © operation but 
does not involve Hermitian conjugation. Unlike in the 
case of a single scalar field, the constraint equations (B39) 
and (B40) are not algebraic and cannot be solved easily. 
The solution of the algebraic part of (B39) and (B40) 
is not a solution of the complete differential equations. 
Therefore, we conclude that even to linear order in the 
gradients the spectrum of the mixing scalar fields de- 
pends on the derivative terms described by the Poisson 
brackets. 



4. Solution in the equilibrium limit 

Since in vacuum and in thermal equilibrium the system 
is homogeneous, isotropic and static the two-point func- 
tions are translationally invariant, i.e. do not depend on 
the center coordinate X. Therefore, in thermal equilib- 
rium and in vacuum, the Poisson brackets (B32) vanish. 
In the following, we discuss the solutions of the constraint 
equations obtained by neglecting the Poisson brackets. 
However, we do not make use of any relations relying 
on periodic equilibrium boundary conditions. Therefore, 
the formal solutions should also be a reasonable approx- 
imation sufficiently close to equilibrium, when gradient 
contributions are small. In this case, (B40) reduces to an 
algebraic matrix equation. Its solution reads 

G R (X,p) = -n n 1 (X,p). (B41) 

From the definition (B34), relation (B28) and the analo- 
gous relations for the retarded and advanced self-energies 
TIa( R ) it then follows that 

G A (X,p) = -Cl A 1 (X,p). (B42) 



From relation (B31) it follows that 

G p (X,p) = i[-G R (X,p) + G A (X,p)} . (B43) 

Using (B41) and (B42) we then find that in equilibrium 
and in vacuum 

G p (X,p)=i[Cl R 1 (X,p)-Cl- A 1 (X,p)] 

= -n R l A) (X,p)U p (X,p)n A l R) (X,p) . (B44) 

Similarly to (B44), the equilibrium (or vacuum) solution 
of the constraint equation for the statistical propagator 
is given by 

G F (x,p) = -n B 1 (x,p)fi F (x,p)n A 1 (x,p) . (B45) 

Note that unlike in (B44) the order of the retarded and 
advanced components is no longer arbitrary. Let us now 
split fl R (A) m to off-diagonal and diagonal matrices, the 

latter being denoted by Q R ( A )- Next we 
define diagonal matrices 

g R(A) {X,p) = -Q-\ A) {X,p) (B46) 

and 

G F{p) (X,p) = -f2 R 1 (X,p)n Hp) (X,p)f2 A 1 (X,p) , (B47) 

where 77f(p) denotes the diagonal components of the self- 
energy matrices n F ( p ). Explicitly 

IF F{p) (X,p) 

gl ^ p) = 'WU^¥TW^^- (B48) 

In the quasiparticle approximation, i.e. in the limit of va- 
nishing self-energy, the diagonal spectral function (B48) 
reverts to (9). From (B46) we also obtain 

^'"-ppJ^feiF- (B49) 

In the quasiparticle approximation it reduces to (37). 
Combining (B44) and (B45) with (B46) and (B47) wc 
can express the full propagators in terms of the diagonal 
ones: 

G F{p) (x,p) = [i + g R (x,p)n R (x,p)Y l 

x [g F{p) (x, P ) - g R (x,p)ff F{p) (x,p)g A (x, P )] 

x [i + ff^X.p^AX^p)]- 1 , (B50) 

where 77 is the off-diagonal part of II. From (18) and 
the definitions of the Wigner-transforms of the statisti- 
cal propagator and spectral function it follows that the 
corresponding Wightman propagators are given by 

G^(X,p) = G F (X,p)±tG p (X,p). (B51) 

Inverting the resulting (two-by-two) matrices and taking 
into account that the products of n R ^ and g R ( A ) arc 
off-diagonal we arrive at (22). 
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5. Boltzmann equation 

To obtain Boltzmann equations we neglect the (con- 
ventional) Poisson brackets on the right-hand side of 
(B35). In the Boltzmann approximation the field's mass 
disappears in the difference of the diagonal components. 
However, the difference of the masses squared, AM^- = 
Mf — Mj, appears in the equations for the off-diagonal 
terms. Moving these terms to the right-hand side we find 

p a V a G F ( p ) =%([Q.rGf( p ) - G F ( p )£Ia] 

-[fl F{p) G A -G R U F{p) ]), (B52) 

where all the functions are evaluated at the same point 
(X,p) of the phase-space. Substituting the equilibrium 
solutions (B44) and (B45) into (B52) we see that its right- 
hand side vanishes indeed. 

If all the matrices in (B52) were diagonal, it would re- 
vert to two independent systems of Boltzmann equations 
for the statistical propagator and spectral function: 

p a V a gf = \ [ QfW; - WpQ* ] , (B53a) 
p a V a gf = . (B53b) 

Equation (B53b) is consistent with the fact, that for a 
single scalar field the solution for the spectral function 
(B48) is valid up to the first order in the gradients. This 
allows us to employ the Kadanoff-Baym ansatz for the 
statistical propagator, Qp = (f^ i + ^)G p l , and rewrite 
the original system of the Boltzmann equations (B53) 
as a Boltzmann equation for the distribution function 
f^ i . In the case of two mixing scalar fields the situa- 
tion is quite different. The solution (B44) is valid only 
in equilibrium. As follows from (B38), out of equilibrium 
the off-diagonal components of the spectral function re- 
ceive corrections linear in the gradients. Strictly speak- 
ing, this means that we can not introduce a generalized 
Kadanoff-Baym ansatz and reduce (B52) to a system of 
two equations for the one-particle distribution functions. 
In general, to analyze the generation of the asymmetry, 
one has to solve the system of the Boltzmann equations 
(B52) for the full spectral function and statistical prop- 
agator. Only when the off-diagonal components are sub- 
dominant (which is the case in the hierarchical and reso- 
nant regime, but not in the maximal resonant limit) one 
can approximate (B52) by (B53). The Boltzmann equa- 
tions for the diagonal propagators (B53) do not contain 
the self-energy CP-violating parameter. Given that the 
contributions proportional to the vertex CP- violating pa- 
rameter also cancel out in the Boltzmann equations for 
the real scalars [2] , this is an expected property. 

Appendix C: Integration over the contour 

In this appendix we calculate two integrals over the 
closed-time-path contour, shown in Fig. 14, which are re- 
quired for the derivation of the self-energy contribution to 



the CP-violating parameter in the Kadanoff-Baym for- 
malism. Let us first consider the integral 

I(x,y)=f@ 4 uA(x,u)B(u,y). (CI) 
c 

Assuming that the two-point functions possess a decom- 
position into statistical and spectral components similar 
to Eq. (B16), we find 

x° 

^F(p) (x, y)=J @ A u A p (x, u)B F{p) (u, y) 
o 

y° 

- J@ 4 uA F{p) (x,u)B p (u,y). (C2) 
o 

These formulas are required to obtain the right-hand side 
of (B20). Replacing A p and B p by A R and Ba we can 
rewrite the right-hand side of (C2) as an integral over 
the whole u axis. From (C2) it also follows, that 

U R (x, y)=J 6{u°) A R (x, u)B R {u, y) . (C3) 

Next we consider the integral 

J{x, y) = ff@ 4 u 3> 4 v A(x, u)B(u, v)C(v, y) . (C4) 
cc 

Since J(x, y) can be represented as an integral of a prod- 
uct of I(x,v) and C(v,y) over the contour, its spectral 
and statistical components are given by (C2) with A re- 
placed by I. Furthermore, using (C3) we find 

Jf (p ) (x, y) = - JJ@ 4 u @ 4 v 9(u°) 0(v°) 

x [A R (x, u)B R (u, v)C F{p) (v, y) 
+A R (x, u)B F{p) (u, v)C A (v, y) 
+A F{p) (x,u)B A (u,v)C A {v,y)}. (C5) 

Building the linear combinations of the spectral and sta- 
tistical components we can easily derive the ^ compo- 
nents from (C5). 

We will now calculate the Wigner-transform of (C5) 
in the Boltzmann approximation. That is, in each of the 
functions in (C5) we neglect the deviation of the corre- 
sponding center coordinate from X = X xy . For instance: 

A(x, u) -> A(X XU , s xu ) -> A(X xy , s xu ) . (C6) 

In this approximation the integration over u and v in- 
duces a simple relation between the momenta q\= qi = 
<7 3 = q. Integration over the relative coordinate s induces 
an additional constraint q = p. Therefore we obtain: 

J F{p) (X,p) = -[A R (X,p)B R (X,p)C Hp) (X,p) 
+A R (X,p)B F{p) (X,p)C A (X,p) 
+A F(p) (X,p)B A (X,p)C A (X,p)}. (C7) 

This completes the calculation of the Winger transform. 
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Appendix D: Renormalization 

In this appendix, we derive the renormalization pre- 
scription (43) employed in Sec. Ill for the derivation of 
the CP-violating parameter in the resonant case. For the 
level of approximation considered there, it is sufficient to 
include pcrturbative one-loop mass and field countcrt- 
erms (see Appendix A) 7 for the real scalar field ipi, 

SC = ^SZijd^d^ - ^SMf^j . (Dl) 

Including the counterterms in the Lagrangian (1) results 
in a modification of the Schwinger-Dyson equation (12), 

G- 1 (x,y)=%- 1 (x,y)-tl ren (x,y), (D2) 

7 The renormalization of the full Kadanoff-Baym equations (B20) 
has been discussed in rcfs. [74, 75], see also [76-80]. 



(using matrix notation; see Appendix B 1), where 
U ren (x,y) = U(x,y) - i(5ZU x + 5M 2 )5° (x,y) . 

Proceeding analogously as in Sec. Ill, we decompose the 
renormalized self-energy into diagonal and off-diagonal 
parts, 

fl ren (x, y) = n ren (x, y) + (7T e ")' (x, y) , 

and define a renormalized "diagonal" propagator by 

[Q rm ]-\x, y) = 9-\x, y) - f[ ren (x, y) . (D3) 

In analogy to the steps leading to Eq. (20), we find for 
the full Wightman functions 



(x, y) = G^ n (x, y) - JJ^u^v 6(u°) 9(v°) [& a (x, u)ff^ (u, v)G r / n (v, y) + G^(x, u)ff A {u, v)g r / n (v, v) 

+G R ,(x,u)ff R (u 7 v)g^ n (v,y)] (D4) 
- J^u 8(u°) [G^(x, u)(SZ'a u + (SM 2 Y)g r / n (u, y) + G R (x, u){5Z'U u + {5M 2 )')g^ n {u, y)] , 

I 



where 5Z' and (5M 2 )' denote the off-diagonal parts of SZ 
and 5M 2 , respectively. Using this, we find that the renor- 
malized version of Eq. (22) can be obtained by replac- 
ing Q -+ G ren , as well as fl' R{A) {X,p) -+ fl' R{A) {X,p) - 

SZ'p 2 + (SM 2 )' . The latter prescription coincides with 
Eq. (43) for i^j. 

The diagonal part of the self-energy enters in the renor- 
malized diagonal propagator g ren - see Eq. (D3). Since 
this equation can be split into independent equations for 
each entry on the diagonal, it is analogous to the case 
of a single real scalar field; see e.g. [19]. The solution of 
the corresponding Wigner-transformcd kinetic equations 
for the retarded and advanced propagator reads (up to 
second-order gradients; see also Eq. (B46)), 

g%%\x, P ) = -\p 2 - m 2 ir™ A f(x,p)}-' , (D5) 

where the renormalized self-energies are determined in 
accordance with Eq. (43) for i = j. Since g r p m {X,p) = 
2ImQ R e ^ (X, p) , the spectral function also contains 
renormalized self-energies. Furthermore, using the 
Kadanoff-Baym ansatz g r ^ n ^ 1 = (J^. + ^)g r p en -' a implies 
that also the statistical propagator, and therefore also 
the Wightman functions, are finite. 

Thus, altogether, we find that the prescription (43) 
renormalizes the diagonal as well as off-diagonal compo- 
nents of the relevant propagators. 



Appendix E: Analysis in the full KB formalism 

In this appendix, we derive a time-evolution equation 
for the "baryon" asymmetry within noncquilibrium field 
theory. It is based on the approximate £7(1 ^-symmetry 
b — >• e la b of the toy-model Lagrangian (1), and does not 
require any further approximations beyond the 2PI trun- 
cation. As we will see, in the Boltzmann limit, the re- 
sulting equation for the asymmetry coincides with those 
obtained from the quantum-corrected Boltzmann equa- 
tions discussed above. 

The Lagrangian (1) can be split into _B-conserving and 
B- violating parts, 

C = Cb + Cab , 

where 

C B = C kin - \m 2 ^ - m 2 bb - ^(W , (El) 
C AB = -^ibb-^ibb. (E2) 
The corresponding Nocther current is given by 
J„(a;) - -2i (b^d^x) - d^b(x)b{x)) . 

Its expectation value can be written as (note that we are 
working in the Heisenberg picture) 

3ii( x ) = ( J v( x )) = 2 * (px* - d y »)D F (x,y)\ x = y , 
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where D F (x,y) = \{b{x)b{y) + b(y)b(x)) is the statistical 
propagator of the complex &-field. 

Because of the presence of Lab, the Noether current 
is in general not conserved. Its divergence reads 

0*%(a;) - \ (d x „ + 9, J [2i (d x » - d y ,)D F (x,y)]\ x=y 
= ip x -d v ]D F (x,y)\ x=v . (E3) 

The expression in the last line can be evaluated by us- 
ing the Kadanoff-Baym equations for the complex scalar 
field. These can be derived analogously to the system 
(B20) of Kadanoff-Baym equations for the real scalar 
fields. They read 

v° 

[D x + m 2 (x)]D F {x, y) = J $> 4 z H F (x, z)D p {z, y) 

o 

x° 

- f@ 4 zZ p (x,z)D F (z,y) (E4a) 
o 

y° 

[D x + m 2 (x)]D p (x, y) = J ® 4 z £„(*, z)D p (z, y) . (E4b) 

x° 

We use the notations of Ref. [2] here. In particular, the 
statistical propagator D F and the spectral function D p 
of the complex field are defined analogously to Eq. (B14). 
The upper equations possess an equivalent representation 
in terms of the antiparticle propagators 

D F (x,y) = D F (y,x) , 

D p (x,y) = -D p (y,x) , (E5) 

and antiparticle self-energies Ej? p (defined analogously). 
The Kadanoff-Baym equations for the antiparticles then 
take the same form as Eqs. (E4), except for the replace- 
ment D -> D and £ -» £. 

The Kadanoff-Baym equations for D F and Dp can 
now be used to evaluate the right-hand side of Eq. (E3), 

9"Mx) =i[a x + m 2 (x)] [D F {x, y) - Dp{x, y)] \ x=y 

x° 

= ij 2> 4 z {Y, F (x, z)D p (z,x) - E p (x, z)D F (z,x) 
o 

- flp(x, z)D p (z, x) + £ p (x, z)D F (z, x)} . (E6) 

By using the retarded and advanced propagators 

D R (x,y) = e(x°-y°)D p (x,y) , 

D A (x,y) = -Q(y° - x°)D p (x,y) , (E7) 

as well as analogous definitions for the antiparticle prop- 
agators, we finally find the following equation of motion 
for the toy-baryon current: 

d^j p {x) = -ij$> 4 ze(z°)x (E8) 
x {Y, F (x, z)Da{z, x) + £,r(x, z)D F (z, x) 
- E_p(x, z)D a (z, x) - E_r(x, z)D F (z, x)} . 



So far, the derivation is quite general and the upper 
equations are formally fulfilled even in the exact the- 
ory. To obtain a closed system we have to express the 
self-energies in terms of the two-point functions. In the 
Kadanoff-Baym formalism the self-energies are obtained 
by functional differentiation of the 2PI effective action. 
Approximations can be obtained by truncating the 2PI 
functional, for example at a certain loop order. This 
leads to so-called ^-derivable approximations, which are 
conserving, i.e. they automatically respect the conserva- 
tion laws of the full theory [33, 81]. This means that, 
in the limit of vanishing _B-violation, the baryon current 
obtained in the Kadanoff-Baym formalism is indeed con- 
served at any loop order of the 2PI functional. 

The structure of the 2PI effective action in two-loop 
approximation can be read off the diagram in Fig. 2 (a): 

iT 2 = - \ g m g*J ^x^y G mn (x, y)D 2 (x, y) 

~ 3 9* m 9nJ S> 4 xS> 4 y G mn {x, y)D 2 (y, x) . (E9) 

Differentiating the effective action with respect to the 
two-point function of the complex scalar field we obtain 

E(x, y) = -gtgjGij (x, y)D(y, x) . (E10) 

The statistical and spectral components of the self-energy 
are given by 

T, F (x,y) = -g*g j [G F (x,y)D F (y,x) 

+ \G i J(x,y)D p (y,x)\, (Ella) 
Z p (x,y) = -g*g j [G l ^{x,y)D p (y,x) 

- G»(x,y)D F (y,x)]. (Ellb) 

The upper expressions, together with the Kadanoff- 
Baym equations (E4), as well as the corresponding equa- 
tions for the toy-Majorana fields (B20) and (13), form 
a closed system of equations. Inserting their solutions 
into Eq. (E8) yields the time-evolution of the toy-baryon 
asymmetry 

n B (t) = ^ f d 3 x j (t,x) . (E12) 

Equation (E8) takes memory and quantum off-shell ef- 
fects into account in a very general way, and is valid even 
very far from thermal equilibrium. Furthermore, due to 
the self-consistency it implicitly resums a large number of 
multiparticle scattering processes. However, this feature 
also makes the (numerical) solution of the full equations 
difficult. In a typical leptogenesis scenario one is inter- 
ested in a situation where the deviations from equilib- 
rium are comparably mild. In this case, it is possible to 
obtain simplified evolution equations for the asymmetry 
from Eq. (E8). In the following, we will show that the re- 
sulting equations are equivalent to those obtained from 
the kinetic evolution equations of the particle distribu- 
tion function, provided both are based on the identical 
truncation of the 2PI effective action. 
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1. Consistency with the kinetic approach 

Usually, within the kinetic approach, the baryon asym- 
metry is calculated by evaluating the difference of particle 
and antiparticle distribution functions, 

n b (t) n- b (t) =f$p [hit, k) - f b (t, k)] , (E13) 

where f b (t, k) and f b (t, k) are obtained from Boltzmann- 
like kinetic equations. Equations (E12,E8) for the toy- 
baryon asymmetry discussed above provide a comple- 
mentary possibility to calculate the asymmetry. In this 
context, the question arises, under which conditions both 
approaches are equivalent. In the following, we will show 
that equivalence holds provided (i) one uses the top-down 
or quantum- corrected kinetic equations, which are ob- 
tained from the gradient expansion of the full Kadanoff- 
Baym equations, (ii) the evolution equation (E8) of the 
baryon current is expanded to the same order in the gra- 
dient expansion, and (iii) both sets of equations are based 
on the same truncation of the 2PI effective action. This 
can be seen as follows: By performing a Wigner transfor- 
mation and a first-order gradient expansion of Eq. (E8), 
and using relations (E5), we obtain on the one hand 

d^(X) = -f*fr[E F (X,p)D p {X,p) 

- H p {X,p)D F {X,p) - {S F , D h } PB 
-{Zh,D F } PB \, (E14) 

where the curly brackets denote the usual Poisson brack- 
ets, and we have used the decomposition D R ^(X,p) = 
D h (X,p)±iD p (X,p). 

On the other hand, the quantum kinetic equation for 
Dp obtained from a first-order gradient expansion of the 
Kadanoff-Baym equation (E4a) reads (see e.g. Ref. [2]) 

{oj(X,p),D F (X,p)} PB = V p (X,p)D F (X,p) 
-X F (X,p)D p (X,p) + {S F , D h } PB , (E15) 

where u> (X, p) = p 2 — m 2 — S/, (X, p) . Inserting this equa- 
tion into Eq. (E14) yields 

^3M) = / $F 2p a V a D F (X,p) . (E16) 
Using D F (X,p) = D F (X, —p), we can write this as 

&*MX) = ftr 2 P a v a [D F (x,p) D F (x,p)] . 

a 

(E17) 

Inserting the Kadanoff-Baym ansatz D F (X,p) = 
[f b {X,p) + \]D p {X,p) yields 

, oo 

d ^M) = f£& Jtr 2 P av * Iff - h\ D P {x, P ), (Eis) 

o 

where we have assumed a symmetric spectrum, D p = 
D p . For a spatially homogeneous system, the frequency 



integration yields the particle distribution function of the 
momentum mode p at time t = X , 

oo „ 

f b (X°,p) ee f% F 2p°f b (X,p)D p (X,p) . 
o 

Thus, using Eq. (E12), we find 

= J a?xd t jo(t,x) 

= Il^d t (f b (t,p)-f b (t,p)) 

= j t (n b (t)-n- b (t)) . (E19) 

Thus, we find that the baryon asymmetry ns{t) inferred 
from the equation of motion of the baryon current co- 
incides with the asymmetry n b (t) — n b (t) obtained from 
the corresponding kinetic approach, under the conditions 
discussed above (this is in accordance with the general 
analysis of Ref. [33]). 



2. Self-energy contribution 

Here, we briefly want to discuss the time-evolution of 
the baryon asymmetry within the 2PI two-loop trunca- 
tion based on the full evolution equation (E8). As dis- 
cussed in Sec. II, the 2PI two-loop approximation cap- 
tures the self-energy contribution to the CP-violating 
decays and inverse decays, whereas the vertex contribu- 
tion requires to take also three-loop diagrams in the 2PI 
functional into account [2]. Thus, the 2PI two-loop case 
provides a possibility to investigate whether an asymme- 
try can be induced by the self-energy-type contributions. 
The latter question has, after some controversial discus- 
sions, been answered e.g. in Ref. [82] within the Boltz- 
mann framework. In the following, we address the same 
question, from a slightly different point of view, within 
nonequilibrium quantum field theory. 

Assuming spatial homogeneity, we can write the total 
baryon asymmetry as 

n B (*) = ^ / d 3 x jo (t, x) = / (0 n B (t, k) . (E20) 

From Eq. (E8), we obtain the following evolution equa- 
tion for the baryon asymmetry in momentum mode k: 

^n B (t,k) = (E21) 
t 

i J dt'{i: F (t,t',\s.)D p (t',t,\s.) - X p (t,t',k)D F (t',t,k) 
o 

- t F (t, t', k)D p (t', t, k) + E,(t, t', k)D F (t', t, k)} , 

where L>(x°,y°,k) = / d 3 xe ik ^- y ^>D(x, y). 

By inserting the 2PI two-loop expressions (Ell) into 
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the upper equation, and using Eq. (B15), we obtain 
|n s (i,k) = 

-2iJdt'{Re(g i9 *)[2G%(D F D p - D F D P ) 

+ G y(D F -D%-l(Dl-Dl))] 

+ ilm(g ig *) [2G% (D F D p + D F D p ) 

+ G y(D F + D 2 F -\(Dl + D 2 p ))]} , (E22) 

where all propagators are evaluated at (f',t, k). In the 
limit of vanishing CP-violation, i.e. for Im(gig*) = 0, 
the Eq. (E22) possesses a baryon-symmetric solution for 
which D F = D F , D p = D p (see Appendix F), and ns ~ 
at all times. For an initial state which is asymmetric, 
the terms proportional to Re(gig*) tend to wash out the 
asymmetry and would, in flat space-time, asymptotically 
drive n B — > 0. 

Here, we are interested in the opposite situation. This 
means, we want to answer the question whether nsif) 
can deviate from zero at some instant of time although it 
vanishes initially, ns(t = 0) = 0. Any change of Ub with 
time requires, of course, nonzero B-violation (i.e. g^ 7^ 0), 
since n B would correspond to a conserved charge oth- 
erwise. Additionally, from Eq. (E22), we see that the 
generation of an asymmetry also requires Im^g*) =/= 0, 
i.e. nonvanishing CP-violation. Finally, a deviation from 
thermal equilibrium is required, since otherwise time- 
translation invariance would require ns{t) to be equal 
to its initial value at all times. Thus, the quantum evo- 
lution equation (E22) obtained from the 2PI two-loop 
approximation is in accordance with the Sakharov con- 
ditions. In general, we therefore expect that ns(t) can 
deviate from n B (0) = for some t > 0, if lm(gig*) 7^ 
and if the system deviates from equilibrium 8 . In order 
to make this statement more quantitative, we consider 
a Taylor expansion of Eq. (E22) around the initial time 
t = 0. The first nonvanishing contribution turns out to 
be of order four, 

—n B (t,k)\ t=0 = 

- 4lm( gi g* 2 )D F d t (G 12 - G 21 )\ t=t , =0 (E23) 

- 2n B (0,k)(Re(g i9 *)G F ' + \ 9l \ 2 D F )\ t=t , =n . 



There are two potential sources for deviations from thermal equi- 
librium, namely the space-time expansion as well as the initial 
conditions. In order to check whether Eq. (E22) can indeed de- 
scribe the generation of an asymmetry, we concentrate on the 
latter source for simplicity. We note that, in order to describe 
equilibrium initial conditions of the full quantum evolution equa- 
tions, it would be necessary to include non-Gaussian correlations 
of the initial state (see e.g. [75]). However, since the latter play 
no role in the kinetic approach (the initial time is formally sent 
to the infinite past when performing a Wigner transformation), 
we do not consider them here. 



Thus, for a symmetric initial state with 7iB(0,k) = 0, 
an asymmetry can arise if Im^!^) 7^ 0, as expected, 
and if d t (G 12 '- G 21 )| t=t , =0 ^ 0. The latter condition is, 
however, not necessary. If d t (G 12 — G 21 )| t=t ' = o = and 
ns(0, k) = 0, the leading term reads 

^n B (t,k)| t = = 4lm( 9l g* 2 )(M 2 - M 2 )D F G 12 \ t=t , =0 . 

Consequently, the full quantum evolution equations sup- 
port the observation that an asymmetry can be generated 
from self-energy-type contributions to the CP-violating 
decays. 



Appendix F: Baryonically symmetric configuration 

The Wightman propagators of the complex scalar field 
are defined by [2] 

D > (x, y) = (b(x)b(y)) = Tr[^ b(x)b(y)} , (Fla) 
D K (x, y) = (b(y)b(x)) = Tr[^ b(y)b(x)\ . (Fib) 

Applying the operator of charge conjugation C to (Fl) 
we obtain 

CD^x^C- 1 =Tr[^ c b(x)b(y)} =D° < (y,x), (F2a) 
CD^x^C- 1 =Tr[& c b(y)b(x)] =D%(y,x), (F2b) 

where the superscript 'c' denotes the charge-conjugated 
quantities. In a baryonically symmetric configuration 
ope _ Therefore in this case D^(x, y) = D^(y,x). In 
turn, for the statistical propagator and spectral function 
this implies 

D F (x,y)=D F (y,x), D p (x, y) = -D p (y, x) . (F3) 

Since the full propagator is related to D F and D p by 

D(x, y) = D F (x, y)-i sign c (x° - y°)D p (x, y) , (F4) 

we conclude that in a baryonically symmetric configura- 
tion D(x, y) = D(y, x) = D(x, y) . 

Differentiating the effective action with respect to the 
two-point function of the real scalar fields we obtain the 
two-loop contribution to the self-energy: 

U^(x,y) = -\g i9 *D 2 {x,y) - \gtg 3 D 2 (x,y) . (F5) 

In the symmetric configuration it takes the form 

U^(x, y) = -\{g l9 ) + g* 9j )D 2 {x, y) , (F6) 

i.e. this matrix is symmetric. This is also true for the 
spectral, statistical, retarded and advanced self-energies 
that can be derived from (F6). 
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Appendix G: One-loop integral 

Here we calculate the one-loop self-energy integrals 
given in Eq. (33), 

L h(p) (X,p) = l6njdR k D F (X, k)D h(p) (X,p- k) , (Gl) 

evaluated for an on-shell momentum p — (Ei,p) where 
Ei = (p 2 + M 2 ) 1 / 2 . We have omitted the superscript 
's' here since the calculation applies also to baryoni- 
cally asymmetric states. We first consider L p . In- 
serting the quasiparticle approximation for the com- 
plex field, see Eq. (7), and the Kadanoff-Baym ansatz 
D F = (\ + f b ) D p , one finds 

L p (X,p) =16tt / dU k [\ + f b (k ■ u)] 2^sgn(fc )(5(fc 2 - m 2 ) 
x 2vsgn(p - k )5((p - k) 2 - m 2 ) , (G2) 

where we have used a manifestly Lorentz-covariant no- 
tation. In the rest-frame of the medium u = (1,0,0,0). 
Note that only the poles with ko > and po — ko > 
contribute since po > Mi > 2m. The integration over ko 
is trivial due to the Dirac-i5 in the first line. In order to 
perform the remaining integrations, we express the inte- 
grand in the rest-frame of the decaying "heavy neutrino" , 
for which p' = (Mi, 0) and u' = (Ei, -p)/M;, 

L p (X,p) = 16^/(0 [§ + fb(k ■ u')] £s5(Mi - 2cj k ) , 

(G3) 

where w fc = (k 2 +m 2 )2 and k-u' = (Lj k E i + 'kp)/M l . The 
Dirac 5- function requires oj k — Mi/2. Using furthermore 
kp = |k||p| cosG, we find 

k-u' = {Ei + r\p\cosO)/2 = E p , (G4) 

where r = (1 — 4m 2 /M 2 ) 1 / 2 . After integrating over |k|, 
one obtains 



L p (X,p)=rJ^[l + 2f b (E p )\ 



(G5) 



If the medium is approximately symmetric, f b ~ f b , we 
can write the result in the form presented in Eq. (34) . 

Let us now turn to Lh- As above, we insert the 
Kadanoff-Baym ansatz for Dp and use the quasiparti- 
cle approximation for D p , yielding 

L h (X,p) = 167T / dU k ^S(k a - Lj k ){ [i + f b (k ■ u)] (G6) 
x D h (X,p - k) + [i + f b (k ■ «)] D h (X,p + k)}. 



In the quasiparticle limit, one can express Dh = ReD_R 
in the form 



D h (X,p) = D h (X,p) 



P 



m 



2 ' 



(G7) 



where V denotes the principal value. As above, we evalu- 
ate the integral in the rest-frame of the decaying particle. 
The result of the calculation reads 



r) 3 uj k M t 
+ [h + h(k-u>)] 



Mi - 2oj k 



V 



Mi + 2oo k 



} 



=Lr(p')+Lr d (x, P ). 



(G8) 



The vacuum contribution, which corresponds to the limit 
fb, f b — > 0, depends only on p 2 due to Lorentz invari- 
ance of the vacuum state. It is logarithmically diver- 
gent. However, after renormalization, only the difference 
L™ c (p 2 ) — L^ ac (p 2 ) appears, which is finite. For the on- 
shell momentum p 2 = M 2 considered here, we obtain 



Lvac 
h 



(M 2 ) 



r In 



frMt 



7^ In 

2ir [i 



V 2m 

2 



- r p \n 



2m 



for m — V , 



(G9) 



where r p = (1 — 4m 2 / fj 2 ) 1 ^ 2 . For an approximately sym- 
metric state, the medium contribution can be written as 

d 3 k f b (E) + f- b (E) V 



Lmed 
h 



(X,P) 



8?r 



(2tt) 2 



M 2 



(G10) 



where E = k-u' 



|p| cos 9] in the limit m — > 0. 



By introducing spherical coordinates and substituting |k| 
for E, we obtain 



jjned 



(X,p) 



dEE(f b (E) + f- b (E))x 
V 



4tt (Ei + |p|cos6) 2 -4E 2 



(Gil) 



After performing the angular integration (and replacing 
Mi -> M™ ed ), we obtain the result stated in Eq. (46). 
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